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ABSTRACT 


This  report  describes  experimental  and  theoretical  investigations  of  the  part¬ 
icle  emission  produced  in  the  laser- surface  interaction.  Preliminary  measure¬ 
ments  on  the  angular  distribution  of  ions  emitted  from  a  tungsten  target  at  power 
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levels  around  14  MW /cm  indicate  anisotropy  of  the  emission,  with  a  strong  en¬ 
hancement  in  directions  near  the  normal  to  the  target.  Ion  and  neutral  molecule 
emission  from  laser- irradiated  sodium  targets  are  described.  Calculations  on 
the  free  expansion  of  the  laser- produced  blowoff  material,  which  include  the 
asymptotic  velocity  of  the  expansion,  yield  pulse  shapes  similar  to  those  observed 
experimentally  if  a  high  initial  temperature  is  assumed.  The  magnitudes  and 
time  histories  of  shock  waves  produced  in  the  target  by  absorption  of  the  laser 
radiation  have  been  calculated. 
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SECTION  I 


INTRODUCTION 

The  investigations  described  in  this  report  are  continuations  of  work  which 

has  been  carried  out  previously.  These  earlier  investigations  have  been 

concerned  with  particle  emission  produced  in  laser-surface  interactions. 

This  work  has  been  described  in  detail  in  previous  reports.  ^  Some 

knowledge  of  the  background  of  this  work,  as  discussed  in  the  previous  reports, 

is  necessary  for  an  understanding  of  the  methods  and  aims  of  the  investigations 

described  in  the  present  report.  A  time-of-flight  spectrometer  has  been  used 

to  measure  ion  emission  produced  by  absorption  of  ruby  laser  radiation  with 
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power  densities  of  the  order  of  50  MW /cm  .  The  ions  were  .round  to  be  mainly 
alkali  metals  with  energies  of  the  order  of  several  hundred  electron  volts.  The 
neutral  molecule  emission  produced  under  similar  conditions  has  been  studied 
with  a  quadrupole  mass  spectrometer  and  found  to  consist  of  thermally  desorbed 
gases  such  as  Hg  ,  CO,  and  C09  along  with  pulses  of  high  energy  neutral 
molecules  having  energies  of  the  order  of  100  eV. 

In  this  report  we  describe  the  extensions  of  these  types  of  measurements.  In 
some  instances  the  same  types  of  measurements  as  described  in  the  previous 
reports  have  been  carried  out  on  a  different  target  surface,  namely  sodium. 

In  other  cases,  the  types  of  measurements  are  new.  We  describe  construction 
of  equipment  for  measurements  of  the  angular  distribution  of  ions  emitted  in 
the  laser-surface  interaction,  and  describe  preliminary  measurements  on  the 
angular  distribution.  The  equipment  construction  and  the  new  measurements 
represent  a  transition  point  in  this  work.  Previously  the  quadrupole  spectro¬ 
meter  and  the  time-of-flight  spectrometer  were  the  basic  instruments  used  to 
investigate  many  different  types  of  surfaces,  and  the  investigations  were 
carried  out  at  relatively  low  laser  powers  per  unit  area.  With  the  introduction 
of  this  new  equipment,  we  shall  extend  the  range  of  measurements  to  higher 
laser  powers  and  also  begin  to  carry  out  additional  types  of  measurements 
designed  to  elucidate  the  phenomena  occurring  in  the  laser -surface  interaction. 


Specifically,  the  experimental  measurements  carried  out  during  this  report 
period  have  involved  measurements  of  the  angular  distribution  of  ions 
emitted  from  tungsten  targets.  It  was  found,  contrary  to  our  previous 
expectations,  that  the  angular  distribution  from  a  tungsten  target  is  anisotropic. 
The  angular  distribution  is  strongly  peaked  in  the  forward  direction  normal 
to  the  target  surface.  It  should  be  emphasized  that  these  results  are  still 
preliminary,  but  it  appears  as  if  the  strongest  emission  is  within  angles 
approximately  15°  from  the  normal  to  the  surface.  The  flux  along  directions 
more  nearly  parallel  to  the  target  surface  drops  off  to  a  relatively  low  level. 

If,  in  subsequent  measurements,  these  preliminarly  results  are  confirmed, 
it  will  have  implications  for  the  previously  calculated  values  of  numbers  of  ions 
emitted  in  the  laser  surface  interaction.  Those  numbers  were  obtained  assuming 
that  the  angular  distribution  was  approximately  isotropic.  We  shall  have  to 
reconsider  these  conclusions  in  tne  light  of  the  present  knowledge  of  the  angular 
distribution. 

Measurements  have  also  been  made  on  a  sodium  target.  Sodium  is  of  interest 
because  the  largest  component  of  the  ionic  emission  from  other  surfaces  was 
found  to  be  sodium  under  our  conditions  of  laser  power.  These  measurements 
were  similar  to  those  described  previously  in  that  they  used  the  time -of -flight 
mass  spectrometer  and  quadrupole  spectrometers.  Results  indicate  large 
amounts  of  sodium  ion  emission  along  with  fragments  due  to  heavy  hydrocarbons 
which  come  from  the  kerosene  in  which  sodium  had  been  stored.  The  neutral 
gas  pulses  also  showed  much  more  hydrocarbon  contribution  than  was  the  case 
with  oth^r  target  materials.  The  high  energy  neutral  molecules  were  also 
observed  with  qualitative  results  similar  to  those  from  other  target  materials 
such  as  carbon. 
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On  the  theoretical  side  some  new  topics  have  been  considered.  These 
include  consideration  of  the  amount  of  space  charge  which  is  produced  in  the 
lase^-surface  interaction,  and  a  reexamination  of  the  expected  ionization 
level  of  the  blowoff  material.  These  considerations,  which  go  into  much 
greater  detail  than  previously,  indicate  that  it  may  be  worthwhile  to  recon¬ 
sider  the  topic  of  absorption  of  laser  light  in  the  blowoff  material  by  the 

process  of  inverse  Bremsstrahlung,  and  the  possibility  of  rapid  heating  of 
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the  blowoff  material  to  high  temperatures.  The  previous  results, 
which  had  indicated  that  this  heating  would  be  low  under  our  experimental 
conditions,  will  have  to  be  reevaluated  in  the  light  of  the  present  considerations 
on  the  space  charge  effects  and  the  amount  of  ionization. 

Another  new  topic  which  is  considered  in  this  report  is  the  calculation  of  the 
stress  profiles  produced  by  thermomechanical  shock  waves  in  the  target 
material.  The  material  expands  as  energy  is  deposited  in  a  thin  surface 
layer  in  a  short  time,  and  a  pressure  pulse  can  be  produced  which  propagates 
into  the  material  at  the  speed  of  sound.  We  have  calculated  the  time  history  of 
the  stress  profile  using  the  PUFF  computer  code.  Results  indicate  that  in 
our  experimental  situation  we  should  obtain  a  sharp  single  pressure  pulse 

g 

propagating  into  the  material  with  an  amplitude  of  the  order  of  6  x  10  dynes/ 

2 

cm  .  The  stress  profile  will  not  be  strongly  dependent  on  the  shape  of  the 
laser  pulse.  Detailed  pressure  profiles  within  the  target  material  are 
presented  for  various  times  after  the  beginning  of  the  laser  pulse. 

We  also  extend  calculations  of  the  adiabatic  free  expansion  of  the  laser- 
produced  blowoff  material,  which  were  described  previously.  In  this  report 
we  introduce  the  asymptotic  value  of  the  velocity  of  expansion  of  a  cloud  of 
material.  We  find  that  we  get  good  agreement  between  the  shapes  of  the 
pulses  calculated  in  this  manner,  and  the  shapes  of  experimentally  determined 
pulses.  A  necessary  condition  to  obtain  sharply  rising  narrow  pulses  of  the 
type  that  we  observed  experimentally  is  that  the  original  temperature  of  the 
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material  before  it  begins  to  expand  is  high,  of  the  order  of  1C  °K.  This 
temperature,  which  produces  adirected  energy  of  expansion  of  the  gas  cloud, 
is  considerably  lower  than  the  temperature  that  would  have  to  be  invoked 
to  explain  the  experimentally  observed  ion  velocities  as  purely  thermal 
velocities.  The  results  of  this  treatment  indicates  that  a  heating  mechanism 
for  the  gas  should  be  investigated  further.  The  picture  that  is  derived  from 
this  result  is  very  suggestive  of  a  heating  of  the  ions  decoupled  from  the 
neutral  gas,  followed  by  a  free  expansion  which  gives  the  observed  directed 
component  of  the  velocities  of  expansion. 

Finally,  we  discuss  a  speculative  model  in  which  the  features  suggested  by 
the  calculations  on  the  adiabatic  free  expansion,  and  the  measurements  on  the 
time -of -flight  spectrometer  and  the  quadrupole  spectrometer;  can  be  synthesized 
to  make  a  complete  picture  of  the  interactions  in  the  laser- produced  blow  off 
material.  This  picture  has  the  advantage  of  providing  consistent  interpretations 
for  the  results  from  the  two  types  of  spectrometers.  It  seems  to  predict  the 
main  features  of  the  observations,  and  suggests  the  heating  mechanism  for  the 
observed  energetic  ions  but  allows  most  <^f  the  material  (e.  g.  ,  desorbed  neu¬ 
tral  surface  gases)  to  remain  cool.  This  picture  must  be  regarded  as  preliminary, 
but  it  does  combine  enough  features  to  be  physically  appealing.  The  main  question 
is  whether  we  can  invoke  a  heating  mechanism  that  will  provide  the  heating  of 
the  ions  to  the  required  initial  high  temperatures  . 

This  report  occurs  at  a  transition  phase  between  the  measurements  conducted 
on  many  different  types  of  materials  at  relatively  low  laser  powers  per  unit 
area,  and  new  types  of  measurements  in  which  the  laser  powers  per  unit  area 
will  be  raised  to  higher  values.  The  results  obtained  previously  on  particle 
emission  and  on  the  phenomena  occurring  in  the  blowoff  material  can  be  used 
to  provide  an  integrated  model  of  all  the  processes  that  take  place.  The 
previous  considerations  have  accounted  for  various  phenomena  essentially  one 
at  a  time.  For  example,  calculations  on  the  expansion  of  the  blowoff  material 
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were  carried  out  separately  from  a  consideration  of  the  heating  occurring 
in  the  blowoff  material.  The  present  report  marks  the  point  at  which  these 
various  types  of  measurements  and  calculations  begin  to  be  synthesized  into 
one  complete  picture  of  the  process  occurring  in  the  laser-surface  interaction. 


SECTION  II 


EXPERIMENTAL  INVESTIGATIONS 
A.  EQUIPMENT  CONSTRUCTION 

Construction  of  the  interaction  chamber  for  measurements  of  angular 
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distribution  of  charged  particles  has  been  completed.  An  Ultek  20 
liter /sec  ion  pump  maintains  the  system  at  a  pressure  of  <  10  Torr. 

The  detectors  used  so  far  have  been  unipolar  detectors  of  the  type  shown  in 
Figure  II - 1.  The  detectors  consist  of  a  stainless  steel  cylinder  closed  at 
one  end  with  a  solid  plate  containing  holes  for  electrical  connections,  and 
a  95%  transparent  tungsten  mesh  at  the  other  end.  Inside  the  cylinder  are 
two  grids  of  the  same  tungsten  weave,  as  well  as  a  solid  stainless  steel  disk 
that  serves  as  the  collector.  One  of  the  grids  is  connected  directly  to  the 
collector  while  the  second  one  is  insulated  from  it  so  that  it  can  be  used 
for  decelerating  potential  measurements.  Each  collector  is  tied  to  an 
independent  BNC  output  connector.  The  decelerating  grids  are  all  connected 
together  and  share  a  common  BNC  connector. 

Previously,  we  had  also  described  construction  of  a  new  larger  ruby  laser 
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to  be  used  in  the  laser-surface  interaction  studies.  Construction  of  this 
item  took  considerably  longer  than  was  anticipated  because  of  a  number  of 
unforeseen  difficulties.  These  included  an  unexpectedly  long  wait  for  the 
roof  prism  with  Brewster  angle  entrance  face.  However,  it  is  in  operation 
now  at  a  power  level  of  approximately  100  megawatts. 

A  partially  assembled  view  of  the  system  is  shown  in  Figure  II -2.  This  is  a 
photograph  of  the  entire  system  with  the  dielectric  coated  elliptical  glass 
reflector  removed.  The  flash  tube  may  be  seen  above  the  ruby  in  its  water 
cooling  jacket.  On  each  side  of  the  cavity  is  a  holder  for  the  cryptocyanine 
Q-switcn.  On  the  left  the  Brewster  angle  roof  prism  is  visible,  and  at  the 
far  right  the  holder  for  the  output  reflector,  which  is  a  sapphire  etalon. 
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This  system  is  pumped  by  a  bank  of  capacitors  with  a  total  capacitance  of 

1200  microfarads,  and  a  maximum  voltage  of  approximately  2700  volts.  The 

total  energy  available  as  input  is  thus  approximately  4400  joules. 

The  maximum  input  energy  was  chosen  to  be  approximately 
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50%  of  the  explosion  energy  of  the  FX-67  flash  tube.  By  using  50%  of  the 
explosion  energy  of  the  flash  tube  as  the  maximum  input  energy,  we  expect 
to  obtain  reasonable  lifetime  from  the  flash  tube. 

The  alignment  of  this  entire  system  is  accomplished  using  a  helium -neon  laser. 
The  angular  orientation  of  the  etalon  is  adjusted  so  that  the  beam  reflected  from 
its  front  surface  is  reflected  directly  back  into  the  helium-neon  laser.  Then 
the  position  of  the  ruby  is  adjusted  so  that  the  beam  passes  directly  down  the 
center  of  the  ruby.  This  may  be  observed  from  the  scattering  of  the  helium- 
neon  laser  beam  within  the  ruby.  The  ruby  is  rotated  so  that  the  reflection  on 
the  Brewster  angle  end  is  a  minimum.  Then  the  cyptocyanine  cells  are 
inserted  in  the  beam  and  they  are  rotated  to  give  a  minimum  reflection  at 
Brewster's  angle.  Finally,  the  roof  prism  at  the  end  of  the  system  is  adjusted 
so  that  the  reflection  from  it  travels  back  into  the  ruby.  This  may  be  checked 
by  inserting  a  pin  hole  in  the  original  helium -neon  laser  beam.  If  the  reflection 
from  the  roof  prism  is  not  exactly  aligned,  the  return  beam  can  be  seen  striking 
the  material  next  to  the  pin  hole.  The  system  is  aligned  until  the  helium - 
neon  laser  beam  traverses  a  re-entrant  path. 

When  this  laser  was  fired  in  the  non-Q-switched  mode,  that  is  with  pure 
methyl  alcohol  in  the  liquid  cells  in  the  cavity,  a  threshold  value  of  about  2000 
volts  on  the  power  supply  was  obtained,  corresponding  to  a  threshold  input 
energy  of  approximately  2,  400  joules.  Thus  there  is  considerable  capability 
for  increasing  the  input  energy  to  overcome  the  losses  in  the  c_vity  that 
result  from  the  bleachable  dye.  The  shape  of  the  pumping  light  pulse  is  roughly 
triangular  with  a  full  width  at  half  height  of  approximately  600  ps.  The  current 
through  the  flash  tube  showed  some  evidence  of  ringing.  The  FX-67  flash  tube 
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is  polarized,  that  is,  current  should  flow  through  it  in  only  one  direction.  If 
current  flows  the  opposite  direction,  damage  to  the  flash  tube  and  shortening 
of  its  life  will  result.  We  measured  the  current  through  the  flash  tube  using 
a  coil  around  the  wire  leading  to  the  cathode.  We  found  that  in  the  original 
stages  of  the  pumping  pulse  the  current  flowed  in  the  correct  direction,  but 
it  went  through  zero  and  there  was  a  reverse  current  approximately  1/3  as  high 
as  the  forward  current.  By  adding  inductance  to  the  system  the  reverse  current 
was  damped  out  considerably.  This  has  reduced, but  not  entirely  eliminated, 
the  reverse  current. 

When  cryptocyanine  is  added  to  the  liquid  cells,  Q -switched  laser  action  is 
obtained  at  an  input  voltage  of  approximately  2400  volts  with  a  single  pulse. 

This  corresponds  to  an  input  energy  of  approximately  3,  400  joules.  The 
pulse  height  can  be  varied  by  increasing  the  cryptocyanine  concentration  and 
working  at  higher  input  voltages.  This  will  be  the  means  of  controlling  the 
amplitude  of  the  output  pulse.  The  pulse  shape  is  shown  in  Figure  II -3,  which 
is  drawn  from  a  photograph  of  an  oscilloscope  trace  of  the  output  of  an  ITT 
FW114  phototube  viewing  a  reflection  of  the  output  pulse.  The  response  time 
of  the  entire  system  (phototube  plus  cable  plus  scope)  is  estimated  to  be  less 
than  10  nanoseconds  so  that  this  should  represent  a  true  picture  of  the  shape 
of  the  output  pulse.  The  base-to-base  width  of  this  pulse  is  approximately 
60  nanoseconds,  and  the  width  at  half  height  is  approximately  26  nanoseconds. 
The  ringing  observable  at  the  end  of  the  pulse  is  probably  due  to  lack  of 
complete  impedance  matching  in  the  termination  of  the  cable. 

The  amplitude  of  these  pulses  has  not  as  yet  been  accurately  measured.  The 
reason  for  this  is  that  at  the  higher  laser  power  our  old  calibration  techniques 
are  not  applicable.  Previously  we  would  fire  the  laser  directly  into  a  TRG 
model  100  ballistic  thermopile  in  order  to  obtain  the  total  energy  in  the  pulse. 
Then  from  the  knowledge  of  the  pulse  shape  we  would  obtain  the  peak  power  in 
the  pulse.  From  these  measurements  the  ITT  phototube  viewing  a  reflection  of 


-10- 


Figure  II- 3  -  Pulse  shape  emitted  from  large  ruby  laser.  Time  goes  left  to 
right  at  100  ns/div.  Drawn  from  photograph  of  scope  trace  for 
clarity  in  reproduction. 


part  of  the  laser  beam  could  be  calibrated  to  give  the  output  of  the  beam  when 
it  was  actually  being  used  for  the  laser-surface  interaction  studies.  However, 
we  fired  the  large  laser  at  a  piece  of  copper  similar  to  that  used  in  the  input 
hole  of  the  TRG  thermopile,  and  found  that  the  copper  surface  was  damaged. 
Thus,  we  could  not  use  the  same  calibration  technique  as  used  previously 
because  this  would  damage  the  thermopile.  Therefore,  new  calibration 
techniques  must  be  devised  to  monitor  the  output  of  this  laser.  Instead  of 
shooting  the  beam  directly  into  the  thermopile,  we  shall  use  the  thermopile 
to  view  the  diffuse  reflection  from  a  magnesium  oxide  block.  In  order  to 
carry  out  this  calibration,  we  must  then  know  accurately  the  fraction  of  light 
reflected  from  the  magnesium  oxide  'dock  into  the  thermopile  ,  This  calibra¬ 
tion  is  not  yet  complete  so  we  cannot  state  exactly  the  power  levels  at  which 
the  new  laser  has  been  operating.  However,  a  preliminary  estimate  indicates 
that  we  have  exceeded  power  levels  of  100  megawatts  with  the  system.  By 
increasing  the  cryptocyanine  concentration  and  the  input  voltage,  we  feel 
that  single  pulses  with  even  higher  power  can  be  obtained. 

When  the  output  of  the  laser  is  focused  with  a  simple,  inexpensive  lens  with 

a  focal  length  of  5  cm,  we  can  produce  a  very  large  spark  in  air  at  atmospheric 

-3  2 

pressure.  The  approximate  focal  area  is  10  cm  .  Since  other  measure¬ 
ments  indicate  that  to  breakdown  air  a  total  power  per  unit  area  of  the  order  of 
11  2 

10  watts/cm  must  be  achieved,  this  measurement  also  indirectly  indicates 

g 

that  the  peak  power  available  is  of  the  order  of  10  watts. 

In  summary  then,  the  larger  new  ruby  laser  construction  has  been  completed, 
and  the  design  goal  of  approximately  100  MW  appears  to  have  been  met 
although  the  measurements  of  the  power  output  are  as  yet  incomplete.  In  the 
next  report  period  we  intend  to  bring  the  laser  into  use  in  extending  the  range 
of  our  studies  of  laser-surface  interactions  to  higher  power  levels. 
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B. 


MEASUREMENTS  OF  ANGULAR  DISTRIBUTION  OF  CHARGED 
PARTICLES 


A  certain  amount  of  preliminary  data  has  been  obtained  on  the  behavior  of  the 
interaction  chamber  as  well  as  on  the  angular  distribution  of  charged  particles. 
Figure  II -4  shows  a  cross-section  of  the  geometrical  and  electrical  arrange¬ 
ments  used  to  obtain  data.  The  output  of  a  Q-switched  ruby  laser  is  aimed  at 
the  target  through  a  5/16"  diameter  hole  cut  in  the  hemisphere  at  30°  from  the 
polar  axis  and  in  a  plane  perpendicular  to  the  plane  of  the  detectors,  ^he 
laser  beam  is  slightly  focused  with  a  simple  lens  to  reduce  its  diameter  and 
prevent  it  from  striking  the  hemisphere.  The  target  is  a  32  mm  diameter 
tungsten  disc  held  1.  6  cm  above  a  grounded  stainless  steel  disc  by  a  ceramic 
insulator.  The  center  of  the  tungsten  disc  is  at  the  center  of  curvature  of  the 
17.  5  cm  diameter  hemisphere.  The  collectors  are  placed  at  15°  intervals 
along  a  meridian.  They  collect  charges  arriving  within  a  5/16"  diameter 
aperture  in  the  hemisphere.  The  solid  angle  they  subtend  at  the  target  is  2 ir  x 

_3 

10  steradians,  which  is  one  thousandth  of  the  solid  angle  subtended  by  the  hemisphere. 

The  data  obtained  so  far  has  resulted  from  irradiation  with  a  1  megawatt  peak 

2 

power  laser  pulse  focused  to  a  spot  of  3mm  diameter  (14  Mw/cm  ).  Two  types 
of  measurements  have  been  made:  In  the  first,  we  compare  the  charge 
collected  by  one  of  the  detectors  with  that  collected  by  the  hemisphere  to  see 
if  the  ratios  of  the  signals  are  equal  to  the  ratios  of  the  solid  angles  of  the 
two  collectors.  In  the  second  set  of  experiments,  we  compare  the  signal  from 
detector  (forward  detector  in  Figure  II - 4 )  to  that  of  collectors  #5  through 
#1.  We  emphasize  that  the  results  presented  below  are  of  a  preliminary  nature 
with  several  features  that  are  not  well  understood;  however,  some  interesting 


-13- 


Figure  II -4  -  Cross-section  of  geometrical  and  electrical  arrangement 
used  for  angular  distribution  measurements 
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qualitative  conclusions  can  be  drawn  from  them. 


Figure  II -5  shows  the  signals  obtained  from  resistors  Rg  and  in  Figure  II -4 
upon  irradiation  of  the  tungsten  target.  In  this  experiment  the  target  was 
grounded  (V=0)  and  so  were  the  decelerating  grids;  R^  and  Rg  were  identical 
50  ohm  resistors,  and  R^.  was  removed  (target  grounded).  Both  traces  were 
triggered  by  the  signal  from  a  919  phototube  monitoring  the  laser  pulse  and 
run  at  5jus/cm.  The  upper  trace  is  the  signal  from  the  hemisphere  (R^)  with  a 
vertical  sensitivity  of  0.  2  V/cm;  the  lower  trace  is  the  signal  from  detector 
#6  (Rg)  at  a  vertical  sensitivity  of  0.  02  V/cm.  The  signal  from  the  hemispnere 
(upper  trace)  can  be  interpreted  as  follows:  During  the  first  microsecond  after  the 
laser  pulse,  the  net  current  is  negative  due  to  the  higher  electron  velocities,  but 
after  the  electrons  are  either  collected  at  the  electrodes  or  recombine,  the 
current  becomes  positive  because  of  the  slower  positive  ions.  Since  no  accelerating 
potentials  are  applied,  there  is  a  large  spread  in  ion  arrival  times,  hence  the 
long  tail  of  the  positive  part  of  the  pulse.  The  signal  from  detector  #6  (lower 
trace)  contains  a  considerable  amount  of  ringing  that  we  have  not  been  able  to 
eliminate  yet;  for  that  reason  the  first  microsecond  of  the  trace  does  not  contain 
useful  information.  Also,  the  negative -going  signal  observed  after  8ps  is 
believed  to  be  ringing.  The  positive  part  of  the  signal  that  starts  at  ^  3.  5ps 
after  the  beginning  of  the  laser  pulse  is  due  to  ions  arriving  at  the  collector  of 
detector  #6;  this  can  be  verified  by  noting  that  the  beginning  of  the  positive¬ 
going  signal  from  detector  #6  (lower  trace)  coincides  in  time  with  the  peak  of 

the  signal  from  the  hemisphere  (upper  trace)  as  predicted  by  a  previous 

(2) 

analysis,  if  the  signals  were  produced  by  the  same  pulse  of  charge  traveling 
from  the  target  to  the  collector,  past  the  hemisphere.  A  further  check  on 
the  interpretation  of  the  data  of  Figure  II - 5  is  the  time  of  flight.  We  know 
from  previous  measurements  on  tungsten  in  the  TOF  spectrometer  that  most 
of  the  emission  consists  of  desorbed  Na+  and  K4  ions  ^  and  if  we  assume 
that  the  large  positive  pulse  from  detector  #6  is  due  to  Na4,  we  arrive  at 
an  ion  energy  of  100  eV,  which  is  in  good  agreement  with  the  energies 
measured  in  the  TOF  spectrometer. 
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Figure  II -5 


-  Signals  from  hemisphere  (top  trace)  and  detector  #6 
(bottom  trace).  Time  base  is  5/is/cm  for  both  traces; 
vertical  sensitivity  0.  2  V/cm  for  top  trace  and  0.  02 
V/cm  for  bottom  trace. 


Integration  of  the  area  under  the  pulses  of  Figure  II -5  yields  a  total  charge 

1 1  *  g 

collection  of  2.  4  x  10  elementary  charges  at  the  hemisphere  and  8.  8  x  10 

from  detector  #6.  Assuming  that  these  preliminary  figures  are  correct, 

we  arrive  at  two  interesting  conclusions:  (a)  The  ratio  of  the  signal  from 

the  hemisphere  to  that  of  detector  #6  is  27,  while  the  ratio  of  the  solid 

angles  subtended  by  those  two  detectors  is  'WOOD.  This  means  that  the 

emission  is  highly  anisotropic  with  a  strong  peak  in  the  forward  direction; 

(b)  The  number  of  charged  pariicles  emitted  per  laser  pulse  is  about 

two  orders  of  magnitude  greater  than  what  we  had  estimated  from  the 

measurements  in  the  TOF  spectrometer.  ^ 

Figures  II-6  a  and  II-6b  compare  the  charge  collection  of  detectors  #5  and  #2 
with  that  of  detector  #6  (see  Figure  II-4)  for  two  different  laser  shots.  The 
vertical  sensitivity  and  time  base  of  all  traces  is  0.  05  V/cm  and  2/.is/cm, 
respectively.  A  potential  of  T0"0  v.  was  applied  to  the  target.  The 
lower  trace  in  both  figures  is  the  signal  from  detector  #6 ;  the  upper  trace  in 
Figure  II-6a  is  the  signal  from  detector  #5,  while  in  Figure  II-6b  it  is  the 
output  of  detector  #2.  The  number  of  elementary  charges  collected  in 
Figure  II-6a  is  2.  6  x  10^  for  the  upper  trace  (#5),  and  2.  8  x  10^  for  the 

9 

lower  trace  (#6).  The  corresponding  numbers  for  Figure  II-6b  are  8.  8  x  10 
10 

and  2.  9  x  10  ,  respectively.  Since  the  signal  from  detector  #6  is  essentially 

constant  in  Figures  II-6a  and  II - 6b,  a  comparison  of  the  signals  from  detectors 
#5  and  #2,  that  are  located  at  15°  and  60°  to  the  polar  axis  respectively,  gives 
a  good  indication  of  the  anisotropy  of  charged  particle  emission.  It  can  be 
seen  from  the  above  numbers  that  the  angular  distribution  of  particles  is 
anisotropic,  peaking  in  the  forward  direction,  but  not  enough  to  explain  the 
anisotropy  suggested  by  the  data  of  Figure  II -5.  However,  before  pursuing 
this  apparent  discrepancy  any  further,  we  must  obtain  further  data  on  the 
angular  distribution  without  applying  a  potential,  to  the  target,  as  well  as  the 
distribution  in  the  plane  perpendicular  to  the  one  where  the  detectors  are 
presently  located.  Such  measurements  will  be  performed  in  the  next  contract 
period. 
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c. 


PROPERTIES  OF  NICKEL  TARGETS 


Because  results  of  our  measurements  on  neutral  gas  emission  from  single 
crystalline  nickel  targets  differed  from  those  obtained  by  David  Lichtman  in 
electron  bombardment  measurements,^  we  made  measurements  on  the 
composition  of  the  nickel  itself  in  order  to  assure  ourselves  that  the  targets 
that  we  were  using  were  in  fact  similar  to  those  used  by  Lichtman  in  his 
measurements. 

We  measured  the  purity  of  the  samples  by  means  of  a  residual  resistivity 
measurement  and  measured  the  single  crystal  orientation  of  the  targets  by 
X-ray  diffraction.  The  purity  measurement  involves  measuring  the 
resistivity  of  the  material  at  room  temperature  and  the  resistivity  at  liquid 
helium  temperature,  and  taking  the  ratio  of  the  two  results.  Theoretically 
for  a  perfect  infinitely  large  lattice  of  nickel  atoms,  the  ratio  of  these 
resistivities  should  be  approximately  62,  000.  Because  of  the  surface 
scattering,  since  the  nickel  targets  are  not  infinitely  large,  and  because  of 
magnetic  scattering  in  magnetic  materials,  this  ratio  should  be  reduced  for 
a  perfect  nickel  lattice  of  finite  extent  to  approximately  2,  000.  If  the  material 
is  not  completely  pure  but  is  of  the  order  of  P9.  9%  pure,  as  claimed  by  the 
manufacturer,  we  would  expect  the  ratio  to  be  of  the  order  of  200. 

The  actual  measured  ratios  both  for  our  specimens  and  for  those  by  Lichtman 
were  in  the  range  15  to  20.  Interpretation  of  this  result  in  terms  of  the  purity 
of  nickel  is  somewhat  dubious,  but  the  relatively  low  value  of  the  ratio 
indicates  that  the  purity  probably  does  not  exceed  99.  9%  .  At  any  rate,  the 
purity  of  the  specimens  that  we  used  was  comparable  to  the  purity  of  Lichtman's 
material. 
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The  X-ray  diffraction  measurements  indicated  that  the  crystals  were  in  fact 
good  single  crystals  since  very  sharp  diffraction  spots  were  obtained  in 
general.  However,  the  orientation  of  the  material  was  not  exactly  as 
specified  by  the  manufacturer.  Some  samples  were  misaligned  by  as  much 
as  8*.  The  (100)  axis  should  have  been  normal  to  the  surface  of  the  targets; 
in  fact,  it  was  tilted  at  angles  up  to  8°  from  the  normal.  However,  the  same 
phenomena  was  observed  in  samples  that  were  used  by  Lichtman,  and  it  is 
felt  that  the  departures  from  perfection  of  the  targets  has  not  caused  the 
differences  between  our  results  and  Lichtman's. 

D.  EMISSION  FROM  SODIUM 

The  ion  emission  from  a  metallic  sodium  target  was  measured  in  the  time -of  - 
flight  spectrometer.  The  interest  in  sodium  arose  because  it  has  been  one 
of  the  largest  components  of  the  ion  emission  from  all  our  previous  targets. 

We  felt,  therefore,  that  by  using  a  sodium  target  we  might  enhance  the  sodium 
peak  to  be  much  larger  than  all  others  and  effectively  provide  a  source  of 
single-mass"  ions. 

The  measurements  were  made  on  reagent  grade  sodium  metal  that  was  stored 

under  kerosene  for  a  period  of  several  years.  A  typical  spectrum  is  shown 

in  Figure  II-7.  The  photograph  is  a  mass  spectrum  obtained  in  the  time-of- 

flight  spectrometer  with  an  accelerating  potential  V  =  +2000  v.  ,  ana  no 

decelerating  potential  (V'  =  0).  The  vertical  sensitivity  and  time  base  of  the 

upper  trace  are  0.  1  V/cm  and  5  /u s/cm  respectively;  the  lower  trace  displays 

the  same  signal  as  the  upper  trace,  but  at  a  sensitivity  of  0.  05  V/cm  and 

time  base  of  2  /u s/cm.  It  can  be  seen  from  the  figure  that  of  all  the  mass 

(1  2) 

peaks  observed  for  m  ?s  40  a.  m.  u.  in  the  emission  from  other  targets,  ‘ 
we  only  have  sodium  left.  However,  several  peaks  corresponding  to  masses 
46,  61,  and  101  a.  m.  u.  are  found  in  abundances  comparable  to  the  sodium 
itself.  Although  we  cannot  make  a  positive  identification,  we  feel  that  those 
peaks  other  than  Na+  are  due  to  heavy  hydrocarbon  fragments  derived  from 
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101  61  46  23 


Figure  II  -  7 


Ion  emission  from  sodium  target.  Time  increases  from 
right  to  left.  Both  traces  are  triggered  by  laser  pulse. 
Vertical  sensitivity  and  time  base  of  upper  trace  are  0.  1 
V/ cm  and  5/us/cm  respectively.  Lower  trace  displays  the 
same  signal  as  upper  trace  but  at  0.  05  V/cm  and  2ps/cm. 
The  TOF  spectrometer  potentials  are  V  =  +2000  v.  and  V' 
=  0  v. 
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mat. 


the  kerosene  in  which  the  metal  was  immersed.  Kerosene  is  a  mixture  of 
hydrocarbons  composed  primarily  of  chains  containing  six  to  ten  carbon 
atoms,  '  hence,  it  is  not  surprising  to  find  such  heavy  fragments.  Also, 
the  sodium  metal  we  have  used  is  very  likely  a  sponge -like  solid  which  would 
readily  absorb  solvent.  The  possibility  of  the  high-mass  peaks  being  due  to 
metallic  contaminants  can  be  ruled  out  on  the  grounds  that  they  do  not  adsorb 
readily  and  their  presence  in  the  bulk  is  specified  by  the  suppliers  of  the 
sodium  metal  to  be  less  than  one  percent.  Hence,  hydrocarbons  are  the  only 
reasonable  source  for  those  peaks. 

In  conclusion,  although  sodium  appears  to  be  a  potentially  useful  material 
for  further  investigations  of  laser -induced  ion  emission,  it  is  going  to  be 
necessary  to  use  material  that  is  not  contaminated  with  organic  solvents. 

The  heavy  contamination  can  probably  be  avoided  by  the  use  of  single  crystal 
material. 

We  have  also  introduced  a  sodium  target  into  the  quadrupole  spectrometer. 

This  was  done  in  order  to  compare  results  from  this  same  target  material 
on  the  two  different  types  of  spectrometer.  We  were  also  interested  in 
determining  whether  we  could  detect  the  base  material,  sodium,  in  measure¬ 
ments  made  on  the  gas  desorption.  Since  sodium  has  a  relatively  high  vapor 
pressure  as  compared  to  other  materials  that  have  been  used  in  the  quadrupole 
spectrometer,  we  felt  that  it  was  possible  that  we  might  observe  the  sodium 
in  the  gas  phase.  The  results  of  the  measurements  on  sodium  in  the  quadrupole 
spectrometer  are  qualitatively  very  similar  to  those  from  the  other  metals 
studied  previously. 

In  measurements  on  the  emission  from  the  gas  phase  we  again  found  large 
peaks  corresponding  to  hydrogen  and  to  carbon  monoxide.  In  addition,  large 
amounts  of  hydrocarbons  were  emitted  from  a  sodium  target.  This  is  to  be 
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expected  since  the  sodium  target  was  immersed  in  kerosene  prior  to  its 
introduction  into  the  vacuum.  Large  mass  peaks  at  masses  14,  15,  16,  26, 

27,  29,  and  30  were  observed  because  of  hydrocarbon  emission.  The  amount 
of  hydrocarbons  desorbed  from  the  sodium  target  is  much  greater  than  the 
amounts  observed  from  other  target  materials.  The  relative  amounts  of  the 
various  species  are  given  in  Table  II - 1.  From  other  materials  such  as 
tungsten,  detectable  hydrocarbo  n  signals  were  observed.  However,  they 
were  always  much  smaller  relative  to  masses  2  and  28.  With  the  sodium 
target,  the  amounts  of  hydrocarbon  are  comparable  to  the  amount  of  mass  28. 
In  fact,  the  emission  of  mass  28  undoubtedly  contains  a  component  of  CgH^ 
as  well  as  CO.  Another  distinctive  feature  is  that  the  largest  single  peak 
here  is  molecular  hydrogen. 

We  macL-  a  careful  search  for  emission  of  sodium  as  neutral  gas,  but  were 
unable  to  find  any  conclusive  evidence  of  a  detectable  amount  of  neutral 
sodium  emission.  However,  even  if  considerable  amounts  of  sodium  gas 
were  produced,  it  would  stick  on  the  first  surface  that  it  struck,  so  that 
the  measured  sodium  gas  pulse  from  a  sodium  target  would  not  be  large. 

The  only  neutral  sodium  which  would  reach  the  detector  would  be  that  which 
follows  the  direct  line -of- sight  path  from  the  target  through  the  ionization 
chamber  and  the  quadrupole  region  and  into  the  detector. 

Because  of  the  presence  of  considerable  amounts  of  oxidation  on  the  surface 
of  the  sodium,  we  attempted  to  heat  the  target.  When  the  target  was  heated 
by  passing  an  electrical  current  through  it,  the  surface  appearance  of  the 
target  changed  and  there  were  changes  in  the  dimensions  of  the  target.  When 
the  sodium  was  originally  placed  in  the  target  holder,  it  was  clean  and  bright 
in  appearance,  having  a  somewhat  silvery  color.  It  was  transferred  directly 
from  the  kerosene  in  which  it  was  stored  into  cyclohexane,  which  was  found 
to  evaporate  more  slowly  from  the  surface  of  the  sodium,  giving  it  protection 
from  oxidation  during  the  time  of  transfer  from  the  original  container  of 
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TABLE  II -1 


RELATIVE  INCREASE  OF  VARIOUS  MASS  SPECIES 
ABOVE  BACKGROUND  -  SODIUM  TARGET 


Mass  Relative  Increase 

Number  Species  (arbitrary  units) 


2 

14 

15 

16 

17 

18 

24 

25 

26 

27 

28 

29 

30 


H2 

CHr 

l 

CH, 

CH^ 

OH 


h2° 


c2h 

C2H2 

C2H3 

C2H4+  CO 

C2H5 

C2H6 


300 

7 

18 

20 

4 

20 

1 

4 

24 

7 

50 

6 

1 
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kerosene  into  the  spectrometer.  While  the  sodium  target  was  being  fastened 

in  the  spectrometer,  the  spectrometer  was  flu  ied  with  helium  gas.  By  the 

time  the  system  was  pumped  down,  the  surface  of  the  sodium  target  had 

assumed  a  light  gray  appearance.  This  was  due  to  oxidation  during  the  time 

of  transfer.  By  the  time  the  system  was  pumped  down  to  a  vacuum  of  about 
-7 

10  Torr,  the  sodium  surface  stopped  changing  appearance  and  remained 
light  gray.  When  the  sodium  target  was  hit  by  the  laser  beam,  the  surface 
film  was  removed  and  the  spot  that  had  been  struck  again  had  a  bright 
silvery  appearance. 

When  the  target  was  electrically  heated,  the  light  gray  surface  film  grew 
thicker  and  darker.  The  spot  that  had  previously  been  hit  by  the  laser  beam 
became  a  light  gray  color.  The  darker  gray  areas  under  this  condition 
required  4  or  5  shots  from  the  laser  beam  before  their  original  silvery 
appearance  could  be  restored.  After  repeated  heating  the  sodium  target 
began  to  flow  and  deform  badly.  It  was  apparent  that  there  were  large 
dimensional  changes  in  the  target  and  it  was  barely  hanging  in  the  target 
holder.  The  surface  layer  had  become  quite  thick,  and  the  surface  color 
had  become  brown  with  a  touch  of  yellow.  The  surface  film  had  the 
appearance  of  a  crust  with  a  scaly  appearance.  The  deformation  and  the 
flow  of  the  target  became  so  complete  that  this  experiment  had  to  be 
terminated. 

When  the  target  was  heated,  there  was  a  detectable  amount  of  sodium  gas 
in  the  system,  for  several  minutes  after  the  electrical  current  v/as  turned  off. 
However,  the  background  pressures  for  other  gases  in  the  system  also  rose, 
and  it  was  not  possible  to  obtain  reliable  data  on  the  laser  emission  after 
the  target  had  been  heated  until  the  system  was  allowed  to  pump  down. 

Referring  back  to  Table  H-l,  we  see  the  relative  effect  of  the  active  chemical 
nature  of  the  sodium  surface  as  compared  to  the  surfaces  of  the  targets 
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previously  studied  in  the  quadrupole  spectrometer.  We  find  relatively  larger 
amounts  of  F^O  on  the  sodium  surface,  for  example.  This  indicates  how 
sodium  chemically  combines  with  water  vapor  from  the  air.  A  relatively  large 
amount  of  hydrogen  is  observed.  This  is  probably  due  to  the  decomposition 
of  t.^e  various  hydrocarbons,  since  we  see  many  hydrocarbon  fragments.  The 
two  groups  CH^  and  C£Hg  with  their  fragments  are  represented  extensively 
in  Table  II- 1.  The  largest  single  component  identified  here  is  mass  26, 

There  are  also  additional  hydrocarbon  groups  around  masses  40  and  50  due  to 
Cg  and  types  of  molecules,  but  we  have  not  investigated  these  in  so  much 
detail  as  the  lighter  hydrocarbon  fragments. 

The  high  energy  neutral  molecule  emission  from  a  sodium  target  was  also 

studied.  This  is  similar  to  that  observed  from  the  carbon  targets  as  described 

previously  \  The  major  portion  of  the  emission  consisted  essentially  of 

a  double  spike  similar  to  the  spikes  as  shown  in  Figure  III-2  of  the  last  semi- 
(4) 

annual  report.  A  typical  example  of  the  high  energy  neutral  molecule 
emission  from  the  sodium  target  is  shown  in  Figure  II-8.  This  shows  clearly 
the  double  spike  feature  as  mentioned  above.  There  is  also  an  earlier  pulse 
which  is  a  feature  not  observed  on  the  high  energy  neutral  molecule  pulses 
from  carbon.  We  have  been  unable  to  identify  the  nature  of  the  first  pulse 
satisfactorily.  The  upper  trace  of  Figure  II - 8  shows  the  same  train  of  pulses 
on  a  slower  time  scale  An  interesting  feature  is  that  in  this  picture  the 
laser  double  pulsed,  and  on  the  upper  trace  we  may  see  the  same  structure 
consisting  of  an  ultraviolet  pulse  simultaneous  with  the  laser  pulse  plus  the 
three  pulses  due  to  high  energy  neutral  molecules,  repeated  again  about  half 
way  through  the  trace. 

Qualitatively,  there  is  nothing  particularly  distinctive  about  the  high  energy 
neutral  molecale  emission  from  sodium  targets.  The  general  features  are 
very  similar  to  what  has  been  described  earlier  for  other  target  material, 
although  there  are  persistent  and  interesting  differences  in  the  details.  By 
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Figure  II  -  8  -  High  energy  neutral  molecule  pulses  from  sodium  target. 

Time  increases  from  left  to  right.  Upper  trace:  20/us/div. 
Lower  trace:  5/us/div.  Laser  power  approximately  66 
megawatts /cm2.  Upper  trace  has  inverted  polarity. 

It  shows  the  effect  of  second  laser  pulse. 
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taking  the  ratios  of  the  arrival  time  of  the  three  neutral  molecule  pulses  as 
shown  in  Figure  II-8,  we  find  that  the  ratios  are  such  that  the  second  and 
third  pulses  are  consistant  with  the  pulses  being  due  to  CO  and  CO2  with 
energies  of  the  order  of  100  eV.  This  is  .he  same  interpretation  as  arrived 
at  for  the  carbon  targets.  However,  the  first  somewhat  lower  broader  pulse 
is  troublesome  in  its  interpretation.  If  this  is  hydrogen,  the  masses 
corresponding  to  the  later  higher  spikes  would  be  uncomfortably  small. 

They  would  fall  in  the  8  to  12  atomic  mass  unit  range.  However,  if  we 
identify  the  second  pulse  as  CO,  then  the  mass  of  the  first  pulse  would  be 
approximately  7,  which  would  be  an  unexpected  value. 

Figure  II -9  shows  data  on  the  arrival  time  of  the  second  pulse  (tentatively 

identified  above  as  CO)  as  a  function  of  laser  power  per  unit.  area.  Here  the 

variation  is  somewhat  less  than  what  we  observed  for  the  case  of  CO  emitted 

from  the  carbon  target,  as  v/e  can  see  by  comparing  the  results  to  those 

(4) 

in  Figure  III-4  of  the  last  s^ mi -annual  report.  The  results  shown  in 
Figure  II- 9  display  a  considerable  amount  of  scatter,  but  there  is  a  definite 
trend  toawrds  decreasing  values  of  the  arrival  time  as  the  laser  power  per 
unit  area  increases.  The  numbers  on  the  ordinate  on  the  right  hand  side  of 
the  scale  refer  to  the  energies  in  electron  volts  if  the  mass  is  28  atomic 
mass  units.  These  numbers  are  derived  from  the  corresponding  times  of 
flight  over  the  distance  from  the  target  to  the  multiplier  detector.  As  the 
laser  power  per  unit  area  increases  the  amplitudes  of  the  high  energy 
neutral  molecule  pulses  are  monotonically  increasing. 
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SECTION  III 


THEORETICAL  INVESTIGATIONS 

A.  ANALYSIS  OF  SPACE  CHARGE  SHEATH  AND  USE  OF  THE  SAHA 

EQUATIONS  FOR  DETERMINING  EXTENT  OF  IONIZATION  IN  A 
PLASMA 

Our  previous  work  has  not  considered  in  any  detail  the  space  charge  which 
would  be  associated  with  the  heated  target  surface.  The  purpose  of  the  analy¬ 
sis  in  this  section  is  to  estimate  the  electron  density  in  and  the  thickness  of 
the  space  charge  sheath. 

Figure  III- 1  shows  a  potential  diagram  appropriate  to  the  case  under  consider¬ 
ation.  In  the  convention  used  here  more  negative  potentials  are  shown  at  higher 
levels  so  that  electrons  tend  to  fall  down  the  potential  contours.  The  diagram 
(not  drawn  to  scale)  shows  the  Fermi  level  as  the  reference  potential  and 
point  A  represents  the  thermionic  work  function.  Potential  difference  B-A 
is  the  space  charge  sheath  of  interest  in  this  section.  If  there  are  significant 
numbers  of  positive  ions  a  plasma  is  formed  with  the  characteristic  potential 
plateau  BC.  The  essentially  charge- free  region  CD  indicates  the  potential 
variation  due  to  the  application  of  a  negative  potential  to  an  electrode  parallel 
to  the  target.  If  there  were  no  plasma  the  profile  would  be  essentially  linear 
from  points  B  to  D. 

The  Richardson  equation  describes  the  outward- directed  electron  emission 
current  density  at  a  given  potential  relative  to  the  Fermi  level 

J  =  AT^  exp  ^  “  e  (  ^  o  ^  ^^tJ 

where  we  have  expressed  the  potential  as  a  sum  of  the  thermionic  work  func¬ 
tion  ill  ,  and  the  potential,  i[i,  in  the  electron  sheath  .j  the  case  where  a 
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NEGATIVE  POTENTIAL 


Figure  III-  1  -  Potential  diagram  for  a  parallel  plate  system  with  emitting 
target  and  a  negative  potential  on  the  second  electrode. 
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strong  ion- accelerating  field  is  applied,  essentially  all  the  electrons  will  be 

returned  to  the  surface  or  remain  in  the  sheath.  Right  at  the  surface  there 

2 

will  be  a  current  d.r.sity  AT  exp{-  e^/kt)  directed  away  from  the  surface 
and  one  of  equal  size  directed  toward  the  surface. 

The  distribution  of  the  x- component  of  the  velocity  of  the  emitted  electrons 

(7) 

is  assumed  to  be  given  by 

dN/NQ  =  (m/2ffkT)1/2exp  (-m  vx2/2kT)  dvx 

At  the  surface,  then,  the  average  x- component  of  velocity  will  be  given  by 

v  =  (1/N  )  fv  dN 
X  O  J  X 

=  (m/2  t  kT)1  l2Jo  V  x  exp(  -  mvx2/2kT)  dvx 
=  (kT/2m7r)1/2 

The  electron  density  at  the  surface  would  then  be 

n  -  2J/ev 
e  x 

=  <87rm/e2k)1/2  AT3/2  exp  {-e  iQ/kT) 

assuming  the  reflected  electrons  have  the  same  velocity  distribution  as  the 
emitted  electrons.  The  electron  density  will  be  a  maximum  at  the  surface 
and  decrease,  roughly  exponentially,  as  the  distance  from  the  surface  in¬ 
creases.  In  the  case  where  there  was  no  plasma  it  would  decrease  to  near 
zero  but  if  there  is  a  plasma  it  would  decrease  to  the  plasma  density.  The 
scale  of  the  sheath  thickness  will  be  determined  by  the  potential  distribution 
in  the  sheath  described  by  Poisson's  equation.  The  extent  to  which  any  electron 
penetrates  the  sheath  depends  On  the  kinetic  energy  of  the  electron.  The  solution 
of  Poisson's  equation  together  with  the  equation  for  ng  would  yield  the  potential 
and  density  distribution  within  the  sheath.  The 


solution  to  this  rather  difficult  problem  is  not  required  in  detail  for  our  pur¬ 
poses  since  the  scale  of  the  sheath  thickness  would  be  indicated  by  the  wel1- 

(8) 

known  Debye  length  ' 

1  /  2 

h  =6.9  (T/n  )  '  ,  in  cgs  units. 

For  tungsten  at  3500°K  the  electron  density  at  the  surface  due  to  thermionic 

14-3  -  5 

emission  is  calculated  to  be  3  x  10  cm  and  the  sheath  thickness  2x10  cm. 

It  is  interesting  to  note  that  as  the  temperature  of  the  surface  increases,  the 

total  number  of  electrons  in  the  sheath  will  increase  but  that  most  of  them  will 

lie  nearer  the  surface. 

The  above  result  for  sheath  thickness  is  considerably  smaller  than  the  assump- 

(4) 

tion  made  previously  that  the  sheath  dimension  would  be  given  by  the  electron 
thermal  velocity  multiplied  by  the  time  since  emission. 

Two  further  points  should  be  made  regarding  the  sheath  formation.  First,  the 
largest  component  of  the  thermionic  work  function  is  due  to  an  electron's  own 

(7) 

image  force  which  it  must  overcome.  As  indicated  by  the  diagram  in  Fig'.re 

III- 2,  the  range  of  this  force  is  about  5x10  cm  and  therefore  the  "surface" 

in  the  analysis  above  is  at  about  this  distance  from  the  physical  surface.  The 

electron  density  within  the  range  of  the  image  forces  is  much  higher.  Second, 

the  above  is  a  steady-state  analysis.  The  time  required  to  reach  the  steady- 

state  would  be  approximately  the  time  for  an  electron  to  travel  the  sheath 

- 12 

dimension.  For  the  example  above  this  would  be  about  3x10  sec. 

The  surface  electron  density,  as  calculated  above,  rises  rapidly  with  temp¬ 
erature,  so  that  high  charge  densities  near  the  surface  can  be  attained;  for 

17  -  3 

example,  ng  reaches  10  cm  at  T  =  5340°K. 

Next  we  will  examine  the  use  of  the  Saha  equations  for  calculating  the  degree 
of  ionization  in  a  partially  ionized  gas  or  the  stage  of  ionization  (singly,  doubly. 
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Figure  III- 2  -  Diagram  showing  the  energy  conditions  near  the  surface  of 
tungsten. 
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triply,  etc.  )  in  a  fully  ionized  gas.  The  Saha  equations  often  predict  results 
which  conflict  with  our  intuition,  thus  causing  some  confusion  and  uncertainty. 
The  purpose  of  this  section  is  to  clarify  what  the  Saha  equations  are  and 
through  a  simple  example  demonstrate  a  rule-of-thumb  which  we  may  use  in 
place  of  our  intuition,  should  it  be  fallible  in  this  case. 


We  are  concerned  here  with  calculating  the  population  density  in  the  various 

quantum  states  of  a  gas  in  thermodynamic  equilibrium.  The  derivations  are 

appropriate  tc  cases  of  complete  thermodynamic  equilibrium  and  fortunately 

to  cases  where  only  local  thermodynamic  equilibrium  (LTE)  prevails.  For 

a  discussion  of  the  conditions  under  which  LTE  prevails  in  laboratory  plasmas, 

(9'. 

reference  can  be  made  to  Griems1  book.  Chapter  6  of  this  book  forms  the 
basis  of  the  following  discussion. 


Figure  III- 3  shows  a  schematic  energy  level  diagram  for  a  certain  chemical 
species  in  ionization  stage  z  -  1  (z  =  1  for  atoms,  z  =  2  for  singly  ionized 
ions,  etc.  ).  The  ratio  of  densities  in  quantum  states  with  energy  levels 
En  and  Em  is  given  by  the  well-known  equation 


Nn_1  =  C1  e*P<-En/kT> 

N*'1  gz~l  exp(-E  /kT) 

m  m  r  m 


(III-  1 ) 


2-  |  1 

where  gn  and  gm  are  the  statistical  weights  (the  number  of  states  at  the 
same  energy  level) of  the  energy  levels  E^  and  E  .  For  most  energy  levels 
the  statistical  weights  are  given  by  the  number  of  orientations  that  the  total 
angular  momentum  can  assume,  i.  e. , 


g 


n 


2J  4-  1 
n 


The  problem  now  is  to  generalize  Equation  (III-  1 )  to  include  the  density  of  the 
free  electrons  associated  with  the  ion.  Referring  again  to  Figure  III- 3,  the 
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population  of  free  electrons  dNg  in  the  energy  interval  dE^  with  kinetic 
energy  is  given  by  an  application  of  Equation  (III-  1 )  as  follows: 

dNe  =  dge  expL_-(En+Ei;+Ek)/kT] 


NZ"  1  g2'1  exp{-  E  / kT) 

n  6n  r  n 


(III- 2) 


^p[-'K  ♦  Ek)/kT] 

Sn 


where  E^  is  the  binding  energy  of  electrons  in  energy  level  E^,  and  dgg  is 
the  statistical  weight  of  the  free  electron  states  in  the  energy  interval  dE^. 
The  derivation  of  dgg  is  given  by  Griem.  The  result,  in  MKS  units,  is 


(III-3  ) 


Note  that  the  result  is  in  terms  of  the  statistical  weight,  g. ,  and  the  density, 
Nj,  of  the  ground  state  of  the  next  higher  ionization  stage.  Let  us  also  note 
here  that  while  standard  quantum-mechanical  procedures  are  used  in  the  der¬ 
ivation  of  Equation  (III- 3),  it  is  this  result  which  leads  to  the  surprising  pre¬ 
dictions  of  the  Saha  equations. 

Substitution  of  (III- 3)  into  (III- 2)  and  integration  easily  leads  to 


NeNl 


mkT 


exp  (-E  n/kT) 


A  more  useful  form  of  the  equation  is  the  one  which  expresses  the  relation¬ 
ship  between  the  total  densities  of  two  successive  ionization  stages  and  the 
free  electron  density. 


N  N 
e 


2Zz(T) 


(T)  l  2n4\ 


mkT 


exp  - 


Ez_1  -  AE2"1 
oo  oo 


(III- 4) 


In  this  equation  the  Z(T)  are  the  partition  functions  and  are  given  by 


Z(T)  =  s  Sn  exP<“En/kT) 


(III- 5) 


In  a  plasma,  the  ionization  energy  is  reduced  by  a  (small)  amount  AE^ 
due  to  the  Coulomb  interaction  of  the  plasma  in  which  the  ion  is  formed.  This 
reduction  in  ionization  potential  limits  the  (otherwise  infinite)  number  of  terms 
in  the  sum  of  Equation  (III- 5)  so  the  partition  functions  may  be  evaluated.  The 
dominant  term  in  the  partition  function  is  the  statistical  weight  of  the  ground 
state  of  the  configuration  in  question. 

In  the  general  case,  for  a  given  chemical  composition  (i.  e.  ,  known  densities 
of  atoms  and  ions  N  of  chemical  species  a)  the  various  atom  and  ion  densities 
in  the  plasma  may  be  determined  as  a  function  of  temperature.  In  addition  to 
the  Saha  equations  (III- 4)  one  needs  only  the  equations  of  conservation  of  mass. 


N  =  Z  N 

3  Z  3 


and  macroscopic  neutrality. 


N  =  T  z  N 
e  z,  a  a 


While  many  situations  require  detailed  quantum-mechanical  information  and  the 
use  of  a  computer  to  obtain  suitable  solutions,  we  may  investigate  a  simple  sit¬ 
uation  to  illustrate  a  point.  Let  us  calculate  the  electron  density  as  a  function 
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of  temperature  for  a  system  in  which  there  is  a  total  density  of  N  =  10 

-  3  a 

cm  of  carbon  atoms  and  ions.  (This  is  approximately  eauivalent  to  a 

pressure  of  3  Torr  at  room  temperature.  )  We  will  consider  the  problem 

in  two  temperature  ranges  -  low  temperature  where  ionization  is  not 

complete  and  high  temperature  where  appreciable  multiple  ionization  is 

achieved. 

In  the  low  temperature  range  of  interest  the  partition  functions  for  the  neutral 

(9) 

carbon  atom  and  the  singly  ionized  ion  are  given  by  Griem  ,  page  560, 

Z°  =  10.  79 
Z1  =  6.07 

These  values  may  also  be  obtained  directly  for  the  case  in  which  the  temperature 
is  not  too  high,  so  that  the  first  few  terms  tend  to  dominate  the  series,  simply  by 
carrying  out  the  summation  in  equation  (III- 5)  using  tabulated  values^10^  for  Jn 
and  for  the  varirus  excited  states.  The  first  ionization  potential  for  carbon 
is  11.  27  volts. 

Expressing  temperatures  in  electron  volts  and  densities  in  cm  Equation  (III- 4) 
becomes: 


4.21  xlO21  T3/2  exp(- 11.  27/T) 
F(T) 


In  the  region  of  interest  4N  /F(T)»1  and  therefore 

a. 


A  plot  of  this  equation  is  shown  as  part  AB  of  the  curve  in  Figure  III- 4.  At 
temperatures  above  about  1  eV  there  is  significant  multiple  ionization.  In  the 
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f  -  1 1 2  z-1 

high-temperature  region  we  calculate  the  temperature  T  where  N 
Z  2  "2  2  +  1 

and  N  are  equal.  (N  and  N  will  be  small  in  comparison. ) 

First  we  see  that 


N  =  (z-1/2)  N 
e  z 

If  we  neglect  the  correction  A  E  z  1  and  take  the  ratio  of  the  partition  functions 
of  successive  stages  of  ionization  to  be  unity  (assumptions  which  will  have  little 
affect  on  our  result)  then  Equation  (III- 4)  becomes 


(z-1/2)  =  7.  5  x  104  (TZ_1/2)  3/2  exp  (-E  l/Tz~lf2) 


17  -3 

for  a  total  carbon  density  of  10  cm  .  The  ionization  potentials  E^  for 

carbon  are  tabulated  in  Table  III- 1  along  with  the  corresponding  solutions  for 
tz-1/2 


TABLE  III- 1 

IONIZATION  OF  CARBON 


z 

Ionization 
Potential 
(volts ) 

Stage  of 
Ionization 
(z-1/2) 

Tz  - 1  /  2 
eV 

1 

11.26 

- 

- 

2 

24.  38 

1.  5 

2.  1 

3 

47.86 

2.  5 

3.  9 

4 

64.48 

3.  5 

5.  1 

5 

391.  99 

4.  5 

27 

6 

489.  84 

5.  5 

31 

The  results  of  this  calculation  are  shown  as  points  labeled  1.5,  2.5,  etc.  in 
Figure  III- 4 „  The  curve  through  these  points  would  approximate  the  electron 


density  curve  in  the  plasma. 
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We  see  that  the  gas  becomes  essentially  fully  ionized  when  the  electron  temp¬ 
erature  is  about  1  eV  and  the  atoms  are  completely  stripped  of  the  orbiting 
electrons  when  the  electron  temperature  approaches  40  eV.  A  good  rule-of- 
thumb  is  that  the  electrons  with  a  given  electron  temperature  (and  LTE  at 
least  prevails)  will  cause  ionization  of  levels  up  to  about  an  order  of  mag¬ 
nitude  higher  in  energy. 

Physically  the  ionization  comes  about  through  two  effects: 

1)  There  are  some  high  energy  electrons  in  the  electron  energy 
distribution  with  enough  energy  to  move  an  electron  from  the 
ground  state  to  a  free  state. 

2)  Most  of  the  ionization  would  occur,  however,  as  the  result  of 
a  multi-step  process  wherein  the  bound  electron  is  moved 
through  successively  higher  excited  levels  until  it  is  finally 
free. 

That  the  equilibrium  populations  are  as  indicated  by  the  equations  and  not  as 
one  might  guess,  however,  is  because  of  the  high  statistical  weight  of  the 
.  free-eiectron  states. 

These  considerations  about  the  ionization  level  and  the  depth  of  the  layer  in 
which  the  electronic  charge  is  concentrated  lead  to  revised  values  for  the 
charge  densities  near  the  surface,  and  hence  to  different  plasma  heating 

(2) 

rates,  than  were  empl  yed  in  our  earlier  consideration  of  plasma  heating. 
The  problem  of  plasma  heating  must  be  reviewed  in  the  light  of  the  preced¬ 
ing  paragraphs.  We  will  see  that  for  the  present  at  least  we  must  leave  open 
the  possibility  that  the  plasma  heating  is  due  to  the  inverse  bremsstrahlung 
process. 

(2) 

The  plasma  heating  calculation  took  the  electron  density  as  a  variable  in¬ 
dependent  of  the  ion  density.  The  indication  that  the  space  charge  sheath  is 
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relatively  small  suggests  that  except  where  consideration  is  specifically  restricted 
to  the  sheath  region,  macroscopic  charge  neutrality  should  be  taken  as  a  boundary 
condition.  This  would  alter  somewhat  the  use  of  the  Saha  equations  in  the  earlier 
calculations. 

The  previous  calculations  also  assumed  thermodynamic  equilibrium  in  the 

electron-ion- neutral  system  although  it  was  pointed  out  that  the  heating  would 

(2) 

go  from  electrons  to  ions  to  neutrals.  It  was  shown  in  an  earlier  report 
that  even  the  transfer  of  energy  from  electrons  to  ions  was  slow  enough  to 
cause  some  concern  regarding  the  equilibrium  assumption.  The  heating  of 
the  neutrals  would  lag  further  behind. 

It  should  be  pointed  out  that  thermodynamic  equilibrium  is  not  necessarily 
required  in  order  to  correctly  use  the  Saha  equations.  Local  thermodynamic 
equilibrium  (LTE)  is  a  sufficient  condition  and  in  laboratory  plasmas  when 
electron- ion  collisions  dominate  the  rate  processes  then  LTE  prevails.  Sit¬ 
uations  can  arise  then  in  which  ions  can  be  multiply  ionized  due  to  hot  electrons 
but  the  ion  themselves  can  remain  relatively  cool. 

In  view  of  this  it  opens  the  possibility  of  eliminating  all  the  un- ionized  blow- 
off  material  from  the  plasma  heating  system  and  considering  separately  tbs 
transfer  of  energy  from  the  electrons  to  the  ions.  The  problem  can  then  be 
divided  into  two  parts:  first  describing  the  manner  in  which  the  plasma  elec¬ 
trons  absorb  energy  from  thp  laser  beam  and  then  how  the  energy  is  transferred 
from  the  electrons  to  the  ions  in  amounts  indicated  by  the  experimental  observa¬ 
tion  of  ion  energy. 

The  model  suggested  by  this  reasoning  is  then  that  the  plasma  in  front  of  the  target 
would  be  formed  by  (predominantly  alkali  metal)  ions  which  were  emitted  from  the 
target  in  the  ionized  state  and  thermionically  generated  electrons.  The  electrons 
would  absorb  energy  from  the  laser  beam  and  subsequently  transfer  some  of  it  to 
the  ions.  Only  preliminary  estimates  have  been  made  using  this  model  at  this 
time.  The  calculations  will  be  refined  further  in  the  next  report  period. 
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B.  CALCULATION  OF  LASER-INDUCED  PRESSURE  PULSES  IN  THE  TARGET 


The  absorption  of  laser  radiation  at  a  surface  can  produce  large  pressure 
waves  in  the  target  material.  This  phenomenon  is  of  considerable  interest. 

One  mechanism  by  which  pressure  pulses  could  be  produced  would  be 
evaporation  of  material  from  the  surface  with  subsequent  heating  of  the 
material  and  recoil  of  the  heated  material  against  the  surface.  A  consider¬ 
able  amount  of  work  on  this  phenomenon  has  been  done  previously.  Some 
work  on  experimental  measurement  of  the  pressure  pulses  obtained  in  this  way 
in  comparison  with  a  model  similar  to  that  of  Pyatt^^  has  also  been  reported. 
This  work  has  mainly  been  carried  out  at  high  laser  power  per  unit  area  where  a 
considerable  amount  of  material  is  removed  from  the  surface. 

There  is  also  another  mechanism  by  which  laser  produced  pressure  pulses  can 
be  generated.  This  does  not  involve  removal  of  any  material  from  the  surface 
and  can  occur  at  considerably  lower  laser  power  per  unit  area.  This  effect  can 
be  termed  a  thermomechanical  effect.  As  t^e  laser  radiation  is  absorbed  in  a 
thin  layer  near  the  surface,  the  internal  energy  of  that  layer  increases.  This 
heated  layer  will  expand  by  thermal  expansion.  If  the  heat  is  absorbed  slowly, 
the  expansion  can  take  place  in  such  a  way  that  adjacent  layers  of  material  can 
relax  and  no  shock  wave  will  be  generated.  In  other  words,  the  adjacent  layers 
of  material  have  time  to  move  out  of  the  way.  An  ordinary  thermal  expansion 
of  the  material  will  occur.  If,  however,  the  thermal  energy  is  deposited  in  a 
thin  layer  in  a  very  short  time,  that  is,  in  a  time  much  shorter  than  the  relaxation 
time  for  the  adjacent  material,  a  pressure  wave  can  be  generated.  The  heated 
material  is  trying  to  expand  before  the  adjacent  material  can  relax  and  allow  it 
to  expand.  Thus,  a  pressure  is  generated  on  the  adjacent  material. 

A  rough  estimate  of  the  time  in  which  the  surrounding  material  can  relax 
can  be  obtained  from  the  following  consideration;  relaxation  time  equals  thick¬ 
ness  of  the  layer  divided  by  the  speed  of  sound  in  the  material,  since  the 
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relaxation  which  would  produce  relief  of  the  pressure  pulse  is  propagated 

by  phonons  which  travel  at  the  velocity  of  sound.  For  laser-heated  surfaces 

-4 

the  an,ount  of  material  involved  is  of  the  order  of  10  cm  thick  if  one 

5 

includes  thermal  conduction.  The  velocity  of  sound  is  of  the  order  of  10 

cm/sec.  Thus,  if  a  significant  amount  of  heating  is  experienced  in  a  time 
-  9 

shorter  than  10  seconds,  there  will  be  insufficient  time  for  the  relaxation 

to  jccur  and  the  heated  material  will  be  trying  to  expand  against  the  surrounding 

layers  Thus,  a  pressure  pulse  will  be  produced.  If  the  scale  of  heating  is 

-9 

long  compared  to  10  seconds,  relief  of  the  pressure  can  be  obtained  during 
the  time  of  the  pulse  and  we  would  expect  no  shock  wave. 

This  effect  is  of  interest  in  an  investigation  of  effects  of  laser  radiation. 
Accordingly,  we  have  carried  out  calculations  on  the  magnitude  and  time 
history  of  pressure  pulses  produced  in  our  experimental  situation.  The 
analytical  tool  used  for  these  calculations  was  the  PUFF  computer  code. 

PUFF  is  a  one-dimensional  Lagrangian  program  which  has  been  used 
previous’y  in  t^e  study  of  X-radiation  effects.  The  program  was  generated 
at.  the  Air  Force  Weapons  Laboratory  in  Albuquerque. 

The  PUFF  code  uses  finite  difference  methods  to  solve  the  hydrodynamic 
equations  using  a  slab  geometry  and  assuming  one-dimensional  flow.  One 
may  follow  a  pressure  pulse  moving  through  a  sample  consisting  of  a  number 
of  separate  materials.  The  number  of  different  materials  may  be  chosen  as  an 
input  variable  and  may  be  an  integer  between  1  and  6.  A  grid  or  mesh  is 
superimposed  on  the  materials  in  the  sample.  The  grid  areas  are  called  zones, 
and  the  grid  lines  are  called  zone  boundaries.  The  mesh  is  numbered  increasing 
from  left  to  the  right.  The  radiation  is  considered  as  being  incident  from  the 
left,  and  the  pressure  pulse  starts  from  the  left  and  progresses  through  the 
mesh  to  the  right.  The  main  PUFF  program  is  used  us  a  control  program  to 
call  various  subroutines  which  carry  out  specific  portions  of  the  calculation. 
There  are  six  subroutines  called  GENRAT.  HYDRO,  REZONE,  EDIT,  EQST,  and 
SSCAL. 
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GENRAT  Is  called  once  at  the  start  cf  a  problem,  and  reads  the  input  data 
and  calculates  the  energy  deposition  in  the  material.  The  hydrodynamic 
calculations  are  carried  out  in  the  subroutine  HYDRO.  HYDRO  calculates 
propagation  of  a  pressure  pulse  through  the  mesh  of  zones  by  a  finite 
difference  method.  Starting  with  the  left  zone  and  with  the  initial  conditions 
at  the  beginning  of  the  problem,  the  zone  velocity,  the  position  of  the  zone 
boundary,  and  the  pressure  in  the  zone  are  calculated.  Then  the  same  quantities 
are  calculated  for  the  next  zone  using  as  boundary  conditions  the  new  pressure  in 
the  preceeding  zones.  These  variable  by  variable,  zone  by  zone,  calculations 
continue  until  all  zones  have  been  advanced  to  a  new  time.  This  is  one  cycle  of 
the  program.  Then  starting  with  the  state  of  the  material  at  the  new  time,  the 
same  variables  zone  by  zone  at  an  incremented  time  are  recalculated.  Thus, 
working  cycle  by  cycle,  a  stress  history  of  the  shock  wave  in  the  material  as 
a  function  of  depth  and  time  is  obtained. 

Since  available  computer  time  limits  a  problem  to  a  few  hundred  zones,  the 
REZONE  subroutine  distributes  the  zones  where  they  are  most  needed  by 
reducing  the  zone  sizes  in  the  main  pressure  pulse  and  increasing  the  zone  sizes 
in  less  active  regions  of  the  mesh.  The  criterion  used  is  that  the  zone  sizes 
are  made  smallest  where  the  pressure  gradient  is  largest.  The  EDIT  sub¬ 
routine  is  used  to  output  a  number  of  key  quantities  at  designated  intervals 
during  the  program.  The  SSCAL  subroutine  calculates  the  incr  ease  in  internal 
energy  in  each  zone  due  to  absorption  of  energy  from  the  laser  beam  during 
the  incremented  time  between  cycles.  The  EQST  subroutine  calculates  the 
equation  of  state  of  the  material  from  which  the  pressure  as  a  function  of  other 
variables  is  obtained.  The  general  flow  of  the  problem  is  shown  in  Figure  III- 5. 

It  is  not  the  purpose  here  to  describe  all  characteristics  of  the  PUFF  computer 

(13) 

program  in  detail.  Details  may  be  found  in  an  AFWL  technical  report. 

However,  it  is  worthwhile  to  devote  a  little  discussion  to  the  pressure  calculations 
carried  out  in  PUFF.  PUFF  is  an  elastic-plastic  program.  In  o+her  words. 
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Figure  III- 5  -  Flow  of  PUFF  Code 


for  pressures  below  the  elastic  limit  it  is  assumed  that  the  material  deforms 
elastically,  and  when  the  pressure  is  released  the  material  returns  to  its 
original  undistorted  state.  When  the  yield  strength  of  the  material  is  exceeded, 
however,  plastic  work  is  done  on  the  material  and  the  material  is  permanently 
deformed.  Doing  the  plastic  work  extracts  energy  from  the  pulse  so  the 
pressure  pulse  is  attenuated  by  passage  through  the  material. 

One  must  be  cautious  to  choose  values  for  the  input  mechanical  properties  of 
the  material  that  reproduce  wave  propagation  experiments.  The  equation  of 
state  that  appears  in  PUFF  is  of  the  form: 

P  =  (Cu  +  Du 2  +  Su3)  {1  -  Gul 2)  +  Go  E 

Here  P  is  the  pressure,  C,  D,  and  S  are  input  parameters  that  must  be 
matched  to  wave  propagation  experiments  and  are  characteristic  of  a  parti¬ 
cular  material,  G  is  the  Gruneisen  ratio  of  the  material,  u  is  th?  compression 
of  the  maierial,  u  -  olpQ-l,  p  is  the  density  of  the  material  at  the  current  time 
and  position,  and  o  Q  is  the  original  uniform  density  of  the  material,  and  E  is 
the  internal  energy.  This  equation  of  state  is  matched  to  Hugoniot  which  is 

(14) 

obtained  from  shock  wave  data.  Such  data  may  be  obtained  in  the  literature 
where  shock  wave  data  have  been  projected  to  the  P~u  plane  and  the  fitting  has 
been  done  to  obtain  values  for  C,  D,  S,  and  G  for  various  metals. 

Modifications  in  the  PUFF  code  were  necessary  to  suit  the  problem  to  our 
particular  application.  The  two  main  modifications  were  the  use  of  experi¬ 
mental  absorption  coefficients  characteristic  of  the  monochromatic  laser 
light,  and  changing  the  input  and  output  statements  in  the  program  to  be 
compatible  with  an  H-1800  computer.  The  PUFF  code  was  originally 
designed  to  handle  a  broad  spectrum  of  X-radiation  as  input.  The  absorption 
coefficients  are  based  on  the  formula: 

. 

iUlO  )i~  A^El 
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where  y/o  is  the  mass  absorption  coefficient,  i  refers  to  the  interval  between 
the  l  and  (i+1)1  absorption  edges  of  the  material,  E  is  the  photon  energy, 
and  A.  and  ET  are  constants  for  the  given  interval  for  the  given  material. 

Between  two  absorption  edges  In  (y/o  )  is  a  linear  function  of  in  E.  This 
treatment  handles  the  absorption  coefficient  very  well  at  X-ray  energies  less 
than  the  K-absorption  edge,  and  it  is  also  capable  of  handling  input  radiation 
with  a  broad  energy  spectrum.  However,  for  cur  particular  problem  it  does 
not  give  a  good  picture  of  the  deposition  of  monochromatic  light  at  optical 
frequencies. 

Accordingly,  we  modified  the  program  by  removing  the  entire  section  of  the 
GENRAT  subroutine  where  the  absorption  coefficient  is  calculated  and 
inserting  statements  to  accept  experimental  values  of  the  absorption  coefficient 
directly  as  input  data.  Choice  of  the  ''ppropriate  values  for  these  input  para¬ 
meters  poses  a  problem  in  that,  during  the  scale  of  time  we  are  interested  in, 
the  input  energy  is  conducted  by  thermal  conduction  from  the  layer  of  material 
in  which  it  is  originally  absorbed  to  a  greater  depth.  In  our  first  attempt  to 
use  the  code,  we  used  the  raw  data  for  optical  absorption  coefficients  at  the  ruby 
laser  wavelength.  Later,  for  reasons  to  be  discussed  below,  we  modified  the 
values  of  the  input  absorption  coefficient  to  use  effective  absorption  coefficients 
that  account  for  conduction  of  heat  in  the  materials,  since  PUFF  makes  no 
provision  for  thermal  conduction.  In  the  case  of  X-radiation,  for  which  PUFF 
was  designed,  this  makes  no  difference  since  the  X-rays  are  absorbed  over  a 
scale  of  depths  such  that  thermal  conduction  is  negligible.  For  our  case  heat 
conduction  becomes  important,  as  we  shall  see  later. 

At  this  point  we  made  a  few  trial  runs  of  the  PUFF  code  using  as  input  for  the 
absorption  coefficient  values  of  the  optical  absorption  coefficient  for  ruby 
laser  light.  These  computer  runs  gave  sharp  pressure  pulses  that  propagated 
from  the  surface  of  the  material  to  the  interior.  The  results  appeared  to  be 
physically  reasonable.  However,  there  was  a  considerable  problem  in  zoning. 
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The  original  zor  i  size  is  introduced  as  an  input  variable.  A  necessary 
feature  to  get  physically  meaningful  results  is  that  the  original  zone 
size  be  small  compared  to  the  depth  over  which  the  radiation  is  absorbed. 

Each  zone  must  absorb  only  a  small  fraction  of  the  radiation  incident  on  its 
front  surface.  Since  the  absorption  coefficients  for  optical  radiation  are  very 
large  compared  to  those  for  X-rays,  the  zones  must  be  much  smaller  for  an 
optical  problem  as  compared  Lc  an  X-ray  problem.  This  means  a  large  number 
of  zones  and,  since  each  cone  has  all  its  variables  calculated  in  each  cycle  of 
the  problem,  this  means  that  the  running  time  for  the  program  becomes 
unreasonably  long. 

The  solution  to  this  problem  comes  from  considering  the  fact  that  heat  is 
actually  conducted  by  ordinary  thermal  conduction  a  distance  that  is  large 
compared  to  the  inverse  of  the  optical  absorption  coefficient  even  in  times  of 

t\ 5) 

the  order  of  nanoseconds.  Thus,  we  can  define  an  effective  absorption 
coefficient  which  is  the  inverse  of  the  depth  to  which  the  energy  is  actually 
distributed  during  the  time  of  the  laser  pulse,  considering  both  the  depth 
over  which  it  is  originally  deposited  and  the  depth  to  which  it  '  conducted 
during  the  pulse  duration  This  absorption  coefficient  can  be  considerably 
smaller  than  the  original  raw  optical  absorption  coefficient. 

With  this  formulation  we  attained  meaningful  results  that  are  summarized  in 
Figure  III-6.  Each  curve  in  the  figure  represents  a  pressure  profile  as  a 
function  of  depth  at  a  particular  time  after  the  start  of  the  laser  pulse.  We 
see  that  for  approximately  1  nanosecond  the  height  of  the  pressure  pulse 
keeps  increasing.  At  that  time  the  position  of  maximum  amplitude  begins  to 
move  toward  the  right  through  the  material  at  approximately  the  speed  of  sound. 
The  amplitude  no  longer  builds  up,  but  instead  gradually  decreases  as  the 
pressure  pulse  moves  to  greater  depths.  As  a  function  of  depth,  we  see  that 
we  get  a  very  sharp  pressure  pulse  propagating  through  the  material.  The 
particular  sample  chosen  here  was  taken  as  70y  thick.  This  particular  choice 
of  thickness  is  not  important  to  the  results  since  the  pressure  pulse  would 
merely  continue  to  propagate  if  the  material  were  thicker,  gradually  being 
attenuated  as  it  went  to  greater  depths  in  the  material. 
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ion  of  depth  in  tungsten  target  struck 
>siting  0.3  cal/cm^  in  30  nanoseconds 
F  code.  Each  profile  represents  a 
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Figure  III-6  -  (Continued)  -  Stress  profiles  as  a  function  of  depth  in  tungsten 
target  struck  by  a  flat  laser  pulse  depositing  0.  3  cal/cm2  in 
30  nanoseconds  as  calculated  by  the  PUFF  code.  Each  profile 
represents  a  different  time  after  the  start  of  the  laser  pulse. 
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Figure  III- 6  -  (Continued)  Stress  profiles  as  a  function  of  depth  in  tungsten 
target  struck  by  a  flat  laser  pulse  depositing  0.  3  cal/cm^  in 
30  nanoseconds  as  calculated  by  the  PUFF  code.  Each  profile 
represents  a  different  time  after  the  start  of  the  laser  pulse. 
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Around  15  nanoseconds  after  the  start  of  the  laser  pulse  the  pressure  wave 

reaches  the  rear  surface  of  the  material.  There  it  is  reflected  as  a  tension 

wave.  If  this  tension  wave  is  large  enough,  it  can  exceed  the  ultimate  tensile 

strength  of  the  material  and  cause  spallation  at  the  back  surface.  The  tension 

wave  moves  back  into  the  material,  propagating  toward  the  left.  At  the  front 

surface  it  would  again  be  reflected  as  a  positive  compress -on  wave  which 

would  move  from  left  to  right.  This  pressure  pulse  would  nounce  back  and 

forth  a  number  of  times  between  the  surfaces  until  it  gradually  damped  out. 

The  laser  energy  per  unit  area  integrated  over  the  pulse  in  this  case  was 

o 

chosen  to  be  0.  3  cal/cm"  in  a  flat  pulse  of  30  nanosecond  duration.  These 
values  are  characteristic  of  pulses  used  in  our  experiments.  The  temper¬ 
ature  at  the  tungsten  surface  would  increase  by  about  1350°C  because  of 
absorption  of  such  a  pulse,  so  that  we  would  not  be  in  a  region  where 

vaporization  is  important.  We  see  that  the  peak  value  of  the  pressure  is  of 
8  0 

the  order  of  6  x  10  dynes /cm~.  By  comparison  the  yield  strength  of  tungsten 

10  2 

is  approximately  1  x  10  dynes /cm  ,  and  the  ultimate  tensile  strength  of 

1 0  2 

tungsten  is  approximately  3  x  10"  dynes /cm  .  Thus  for  our  case  we  should 
have  only  an  elastic  deformation  and  we  would  expect  no  spallation  to  occur 
at  the  back  surface  of  the  material.  As  the  laser  power  per  unit  area  is 
increased,  we  expect  this  conclusion  would  be  changed.  First,  as  the  maxi¬ 
mum  stress  became  larger  than  the  yield  strength,  permanent  deformation 
would  occur,  and  as  the  laser  power  per  unit  area  increased  still  further,  we 
would  expect  the  possibility  of  scabbing  of  the  back  surface  of  the  thin  slab 
as  the  ultimate  strength  of  the  material  was  exceeded. 

The  dominant  feature  is  the  large  single  stress  wave  moving  to  the  right 
through  the  material.  The  speed  of  sound  is  calculated  to  be  approximately 
.4  x  10  cm/sec.  The  duration  of  the  laser  pulse  is  30  nanoseconds  but  the 
maximum  amplitude  of  the  stress  wave  stops  increasing  relatively  early  in 
the  laser  pulse.  In  fact,  the  maximum  amplitude  is  reached  around  1.0 
nanoseconds  after  the  beginning  of  the  laser  pulse.  At  this  time  the  maximum 
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stress  is  5.  87  x  10  dynes /cm  .  The  position  of  the  maximum  stress  is  at 
approximately  .  0005  cm  from  the  left  edge.  Previous  to  this  time  the 
maximum  amplitude  of  the  stress  has  been  increasing  steadily  as  the  stress 
wave  moves  to  the  right.  After  this  time  the  amplitude  of  the  stress  wave 
increases  only  slightly  and  then  goes  into  a  gradual  decline.  The  attenuation 
as  the  stress  wave  passes  through  the  material  is  to  be  exper  ed.  Since 
the  maximum  stress  di.es  not  exceed  the  yield  strength  of  the  material, 
no  plastic  work  is  done,  and  we  would  expect  the  wave  to  be  attenuated  only 
slowly  by  the  gradual  expansion  of  the  entire  slab  of  material.  This  in  fact 
is  observed.  Around  15  nanoseconds  the  stress  wave  reaches  the  rear  sur¬ 
face  and  is  reflected  with  about  a  20%  decrease  in  amplitude.  Thereafter  it 
propagates  to  the  left  as  a  rarefaction  wave. 

A  surprising  feature  of  this  whole  process  is  that  there  is  no  succeeding  stress 
wave  produced  by  absorption  of  energy  at  times  greater  than  1  nanosecond. 

Very  early  in  the  laser  pulse  the  maximum  stress  is  achieved  and  thereafter 
the  left  surface  of  the  material  remains  approximately  at  zero  stress  until 
the  rarefaction  wave  returns.  There  is,  of  course,  some  question  as  to 
whether  this  might  possibly  be  an  artifact  of  the  hydrodynamic  program.  We 
have  considered  this  question  carefully  and  have  come  to  the  conclusion  that 
this  represents  the  true  behavior,  at  least  the  true  behavior  of  the  hydrodynamic 
model  used  in  PUFF.  The  reason  that  the  stress  wave  stops  increasing  in 
amplitude  after  approximately  1  nanosecond  is  that  the  position  of  the  maximum 
amplitude  has  moved  out  of  the  region  where  the  energy  is  being  deposited, 
at  least  for  a  flat  pulse.  This  would  probably  be  different  if  a  triangular  pulse 
shape  were  used.  At  the  depth  of  .  0005  cm,  where  the  amplitude  stops  increas¬ 
ing,  the  energy  deposition  integrated  over  the  entire  laser  pulse  is  approxi- 
mately  .45  x  10  ergs/gram,  or  about,  one  quarter  of  what  it  was  at  the  surface. 
As  the  stress  wave  propagates  to  deeper  depths  it  does  indeed  increase  very 
slightly,  but  it  is  passing  through  a  region  where  the  rate  of  energy  deposition 
has  become  very  low.  On  the  other  hand,  to  the  left  of  the  position  where 
the  maximum  amplitude  occurs,  we  can  observe  that  the  hydrostatic  pressure 
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does  indeed  continue  to  increase  with  time  as  further  laser  energy  is 
deposited  in  these  regions.  The  deposition  of  energy  leads  to  an  increase 
of  the  hydrostatic  pressure  through  the  equation  of  state  subroutine.  How¬ 
ever,  from  a  detailed  examination  of  the  profiles  we  deduce  that  the  stress 
to  the  left  of  the  position  of  the  maximum  is  near  zero  even  though  the  hydro¬ 
static  pressure  is  large.  The  reason  for  this  is  that  the  hydrostatic  pressure 
is  almost  exactly  balanced  off  by  a  negative  deviatoric  stress.  The  deviatoric 
stress  measures  the  stress  produced  by  the  expansion  of  the  material.  The 
deviatoric  stress  is  calculated  in  the  HYDRO  subroutine  and  is  a  function  of  the 
expansion  of  the  material  and  the  velocity  with  which  it  is  moving.  The 
deviatoric  stress  is  a  relief  of  the  hydrostatic  pressure  by  motion  of  the 
material. 

Thus  a  physical  picture  of  the  processes  occurring  is  that  at  the  very  beginning 
of  the  laser  pulse  when  the  material  is  at  rest  and  in  its  uncompressed  state  the 
deposition  of  energy  begins  to  produce  a  stress  which  is  relieved  by  expansion 
at  the  left  edge  of  the  material.  Later  when  the  material  is  already  expanding 
fairly  rapidly,  the  expansion  relieves  the  stress  produced  by  the  increment  of 
the  deposited  energy  rapidly  enough  so  that  there  is  no  net  stress  in  the 
material.  It  must  be  remembered  that  for  these  runs  the  total  stresses  are 
relatively  low  compared  to  the  ultimate  strength  of  the  material.  We  nowhere 
exceed  the  yield  strength  of  the  material  which  is  taken  to  be  approximately 
1  x  10  dynes /cm  .  This  situation  is  different  when  larger  amounts  of 
energy  are  deposited  and  where  the  maximum  stress  will  begin  to  approach 
the  yield  strength  of  the  material.  In  that  case  the  stress  can  no  longer  be 
relieved  rapidly  enough  near  the  left  edge  of  the  material. 

An  interesting  feature  of  this  result  is  that  the  stress  profile  is  relatively 
independent  of  the  laser  pulse  after  the  first  nanosecond.  Thus,  we  should 
expect  at  the  rear  surface  shock  waves  of  approximately  the  same  amplitude 
and  duration  even  if  we  were  to  cut  off  the  end  of  the  laser  pulse.  That  is,  a 
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laser  pulse  of  one  nanosecond  duration  should  produce  the  same  shock 
wave  at  the  rear  surface  as  a  laser  pulse  of  30  nanosecond  duration, 
provided  the  power  per  unit  area  delivered  by  the  two  pulses  is  the  same, 
and  the  peak  stress  is  less  than  the  yield  strength.  This  prediction  could  be 
investigated  experimentally  by  using  high  frequency  pressure  transducers 
in  contact  with  the  rear  surface  of  the  sample.  This  also  means  that  the 
artifact  of  introducing  details  in  the  thermal  conduction  process  does 
not  have  a  great  effect  on  the  stress  wave  profile,  as  we  shall  see  below. 

This  occurs  because  the  stress  wave  propagates  out  of  the  region  where 
enei  gy  is  being  conducted  quite  early  in  the  pulse  and  the  thermal  conduction 
will  occur  only  in  the  region  where  the  stress  wave  has  been  relieved 
and  where  the  relief  processes  are  keeping  up  with  the  increase  in 
hydrostatic  pressure  during  the  subsequent  portion  of  the  pulse. 

An  alternative  way  of  lookir  _  at  this  is  to  consider  the  time  for  a  stress 
relief  wave  to  be  propagated  across  the  layer  where  energy  is  being 
deposited.  Stress  relief  propagates  at  the  speed  of  sound  in  the  material, 
leading  to  a  characteristic  time  of  around  one  nanosecond,  as  we  saw 
previously.  This  is  in  fact  the  length  of  time  after  which  the  stress  stops 
increasing  in  amplitude.  Since  the  stress  in  the  layer  in  which  the  energy 
was  deposited  is  relieved  by  expansion  of  that  layer,  and  therefore 
a  pressure  is  exerted  against  the  next  layer,  we  expect  the  stress  to 
propagate  into  the  material,  rather  than  dropping  to  zero. 

Some  light  has  also  been  shed  on  this  process  by  carrying  out  PUFF  code 
runs  for  larger  incident  powers  per  unit  area.  If  the  power  per  unit  area 
is  increased,  the  increase  in  the  peak  amplitude  of  the  stress  wave  increases 
approximately  linearly  with  the  increase  in  laser  power  per  unit  area,  up 
to  the  point  where  the  maximum  stress  amplitude  is  equal  to  the  yield 
strength  of  the  material.  The  shape  and  time  history  of  the  pressure 
pulse  remain  unchanged  up  to  this  point.  As  the  yield  strength  of  the 
material  is  exceeded,  the  maximum  pressure  in  the  sharp  pulse  continues 
to  increase,  but  the  time  history  changes  markedly.  The  first  sharp  pulse 
is  followed  by  other  broader  pulses  of  larger  amplitude.  The  first  sharp 
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pulse  still  propagates  from  left  to  right  in  the  material  as  previously, 
but  now  the  amplitude  dnmps  considerably  faster  at  larger  depths. 

We  interpret  this  as  being  due  to  extraction  of  energy  from  the  pulse 
to  produce  plastic  work  on  the  material.  The  later  succeeding  stress 
waves  that  also  propagate  into  the  material  are  rather  complicated  and 
we  interpret  these  as  caused  by  the  material  near  the  surface  being 
strained  past  its  elastic  yield  point,  and  thus  no  longer  relaxing  fast 
enough  to  relieve  this  pressure  pulse.  Thus,  as  the  laser  energy  is 
raised  enough  so  that  the  elastic  limit  of  the  material  is  exceeded  by 
the  stress  wave,  the  time  history  becomes  considerably  more  complicated. 

We  have  not  pushed  the  investigation  of  the  details  of  these  phenomena 
to  the  point  where  we  feel  we  can  completely  understand  them. 

An  alternative  way  of  presenting  the  results  of  the  PUFF  code  program 

is  shown  in  Figure  III- 7  where  we  present  the  stress  as  a  function  of 

time  at  the  rear  surface  of  a  70 p  thick  specimen  of  tungsten  which  is 

2 

struck  by  a  laser  power  of  0.  3  cal/cm  at  the  front  surface.  For  approximately 
8  nanoseconds  there  is  no  buildup  of  pressure;  then  a  positive  pressure  is 
experienced  as  the  compression  wave  reaches  the  rear  surface.  Around  15 
nanoseconds  after  the  start  of  the  laser  pulse  a  sharp  tension  pulse  occurs 
at  the  back  surface  as  the  original  compression  pulse  is  reflected  as  a 
tension  wave.  This  is  the  time  at  which  scabbing  of  the  back  surface 
would  occur  if  +he  tensile  stress  were  to  become  high  enough.  After 
this  time  the  stress  drops  to  a  relatively  low  value.  At  a  later  time  when 
the  tension  wave  is  reflected  from  the  front  surface  and  returns  to  the 
back  surface,  with  a  relatively  short  period  of  high  negative  stress. 

It  should  be  noted  that  PUFF  is  a  one-dimensional  computer  code. 

The  conditions  of  applicability  here  are  that  the  transverse  dimensions 
of  the  laser  spot  should  be  much  larger  than  the  thickness  of  the  foil  we 
are  considering;  otherwise  the  one -dimensional  approximation  will 
break  down. 
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The  buildup  of  the  compress ional  stress  wave  in  its  early  stages,  at  a  time 
short  compared  to  one  nanosecond,  may  be  seen  from  Figure  III-6.  This 
early  buildup  occurs  at  the  depth  in  which  the  laser  radiation  is  being 
absorbed.  Later,  as  the  peak  pressure  pulse  propagates  to  the  right,  out 
of  the  depth  in  which  the  laser  energy  in  being  deposited,  this  buildup 
stops.  We  get  the  somewhat  surprising  result  that  the  time  history  of  the 
stress  pulse  is  relatively  independent  of  the  power  per  unit  area  in  the 
laser  pulse,  so  long  as  the  peak  stress  does  not  exceed  the  yield  stress 
of  the  material.  The  general  shape  of  the  pressure  pulse  at  the  b;  c  c  surface 
will  be  approximately  the  same  as  that  shown  in  Figure  III- 7,  independent 
of  the  power  per  unit  area  delivered  by  the  laser  pulse.  The  amplitudes  of 
course  will  change  approximately  linearly  with  the  deposited  energy. 

As  an  additional  check  on  the  physical  reasonableness  of  this  interpretation, 
we  changed  the  REZONE  routine  to  allow  the  zones  near  the  surface  always 
to  remain  small.  It  was  thought  possible  that  as  the  peak  pressure  pulse 
propagated  outside  of  the  depth  in  which  the  laser  energy  was  being  de¬ 
posited,  and  the  zones  near  the  surface  became  large,  this  might  affect  the 
physical  nature  of  the  problem.  Therefore,  we  modified  the  REZONE  sub¬ 
routine  so  that  the  zones  in  which  the  laser  energy  were  being  deposited 
would  always  remain  small.  We  found  that  this  made  no  difference  in  the 
pressure  pulses  calculated.  We  interpret  this  as  an  additional  check  on  the 
model  that  we  discussed  above,  and  as  confirmation  that  the  single  sharp 
pressure  pulse  propagating  into  the  material  is  a  representation  of  physical 
reality  rather  than  an  artifact  introduced  in  the  computer  program. 

We  have  also  considered  the  possibility  of  accounting  for  finer  details  of 
the  thermal  conduction  of  heat  into  the  specimen  during  the  laser  pulse.  In  the 
work  discussed  above  we  took  one  single  absorption  coefficient  characteris¬ 
tic  of  the  radiation.  This  in  effect  integrated  the  heat  deposition  over  the  entire 
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laser  pulse.  The  PUFF  code  has  the  option  of  allowing  several  spectra. 

If  we  put  in  a  number  of  different  spectra  to  represent  the  laser  pulse, 
and  have  the  different  spectra  active  at  different  times,  we  can  reproduce 
the  effect  of  thermal  conduction  into  the  material  during  the  laser  pulse 
more  exactly  than  by  a  single  absorption  coefficient  representing  the 
distribution  of  energy  by  thermal  conduction  as  integrated  over  the  entire 
laser  pulse.  For  example,  we  break  the  laser  pulse  up  into  10  intervals  of 
3  nanoseconds  apiece,  representing  a  total  laser  pulse  30  nanoseconds  long. 

In  the  first  three  nanoseconds  we  choose  a  large  absorption  coefficient  to 
represent  how  the  energy  is  being  deposited  near  the  surface.  For  the 
second  interval  we  choose  a  smaller  absorption  coefficient  to  represent  the 
fact  that  heat  has  propagated  farther  into  the  material,  and  so  on,  reducing 
the  effective  absorption  coefficient  for  each  successive  three  nanosecond 
interval.  In  this  way  we  can  simulate  conduction  of  heat  into  the  material 
more  exactly  than  was  done  in  the  procedure  described  above.  The  appro¬ 
priate  modifications  in  the  computer  program  were  made  and  the  calculation 

carried  out  using  effective  absorption  coefficients  for  the  different  spectra 

(15) 

to  represent  results  such  as  have  been  calculated  previously  for  the 
depth  to  which  energy  is  conducted  throughout  the  time  of  the  laser  pulse. 
These  results  were  not  significantly  different  from  the  results  shown  in 
B'igure  III  -  6.  The  peak  pressure,  the  time  history,  and  the  attenuation  of  the 
sharp  pressure  pulse  as  it  moves  into  the  material  wero  approximately  the 
same.  This  means  that  our  results  are  not  critically  dependent  on  the  exact 
form  of  the  absorption  model  chosen.  This  is  reasonable  and  gives  an  added 
confirmation  that  our  results  represent  reality  instead  of  being  a  computer 
artifact. 

C.  EXPANSION  OF  BLOWOFF  MATERIAL 

(4) 

In  the  last  semi-annual  report  we  examined  the  expansion  of  the  laser 
produced  blow  off  material  using  a  model  of  a  one -dimensional  adiabatic  free 
expansion  and  found  that  the  directed  energies  of  the  ions  produced  could  be 

-  61" 


'-  A3*  SfcVrS&i 


higher  than  the  thermal  energies  of  ions  in  an  original  cloud  o’  material 
at  a  characteristic  temperature  Tq  .  After  the  cloud  of  laser  desorbed 
material  is  produced  in  front  of  the  target,  it  will  expand.  If  the  material 
is  hot,  the  original  thermal  energy  will  be  transformed  into  energy  of  expansion. 
By  fairly  simple  equations  we  determined  the  expansion  of  the  gas  and  its  cool¬ 
ing.  The  thermal  energy  appears  at  a  distance  from  the  interaction  as  a 
directed  energy  of  <  xpansion.  '{lie  velocity  of  the  expanding  edge  of  the  gas 
cloud  would  yield  an  upper  limit  to  the  energy  that  would  be  deduced  from 
measurements  on  the  time  of  flight  of  the  particles.  This  energy  is  greater 
than  the  thermal  energy  originally  characteristic  of  the  particle  in  the  gas 

cloud.  Other  workers  in  the  field  have  used  similar  models  for  analysis  of 

(16) 

expansion  of  laser -produced  blowoff  material.  ' 

(4) 

In  the  last  semi-annual  report  equations  were  derived  for  the  motion  of  the 

expanding  edge  of  the  gas  cloud  as  a  function  of  time  and  of  original  temperature 

of  the  gas.  From  these  values  the  gas  temperature  and  the  particle  energy 

corresponding  to  a  particle  traveling  at  the  edge  of  the  gas  cloud  were  also 

derived  as  functions  of  time.  This  was  carried  out  in  a  one -dimensional  model 

for  two  postulated  cases  of  gas  density  profiles.  First  a  uniform  gas  density 

behind  the  edge  of  the  cloud  was  assumed;  second,  a  density  gradient  model 

was  postulated  in  which  the  gas  density  always  remained  highest  at  the  position 

of  the  surface  and  fell  off  to  zero  at  the  position  of  the  expanding  edge.  The 

equations  were  evaluated  numerically  for  the  density  gradient  model,  and  data 

were  presented  on  calculated  particle  fluxes  as  a  function  of  time  at  a  position 

(4) 

1  cm  from  the  surface  for  various  assumed  initial  temperatures  Tq.  It 
was  found  that  for  high  temperatures,  of  the  order  of  100,  000°K  initially, 
steeply  rising  pulses  with  half  widths  of  the  order  of  1  or  2/us  were  predicted. 
These  pulse  shapes  are  very  similar  to  those  actually  observed  from  ions 
in  the  time-of-flight  spectrometer.  If  this  model  has  validity,  it  would  imply 
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that  the  blowoff  material  is  greatly  heated  in  a  short  time  before  it  has  a 
chance  to  expand.  Thus  this  treatment  has  the  built-in  difficulty  that  we 
must  provide  a  heating  mechanism  for  the  gas.  However,  the  directed 
particle  energies  that  are  obtained  are  approximately  seven  times  higher 
than  the  average  value  of  the  one -dimensional  kinetic  energy  associated 
with  the  original  assumed  temperature  Tq  .  This  means,  for  example,  that 
to  explain  the  directed  energies  of  the  order  of  200  eV/ion  one  need  heat  the 
gas  only  to  approximately  30  eV.  Moreover,  this  figure  of  30  eV  is  of  the 
same  order  of  magnitude  as  the  spread  of  energies  in  the  spectrum  deduced 
from  the  pi^se  shapes  observed  on  the  time-of-flight  spectrometer. 

In  the  present  work  we  have  carried  out  similar  calculations  for  the  uniform 
density  model  and  have  also  expanded  the  whole  model  to  take  into  account 
a  more  correct  formulation  of  the  boundary  conditions,  i.  e.  ,  we  have  used 
the  proper  asymptotic  velocity  characteristic  of  an  adiabatic  free  expansion. 


We  will  consider  first  the  case  in  which  we  assumed  that  the  blowoff  material 

had  a  uniform  density  from  the  target  surface  out  to  the  expanding  edge  of  the 

material.  The  treatment  for  this  case  was  developed  in  the  last  semi-annual 
(4) 

report.  It  relies  essentially  on  solution  of  the  equation: 


d 


0 


(III  -  6) 


where  X  is  the  distance  of  the  edge  of  the  material  from  the  surface.  This 
equation  was  derived  neglecting  possible  absorption  of  laser  light  in  the  blowoff 
material.  It  assumed  only  the  existence  of  a  cloud  of  gas  at  an  original 
temperature  Tq.  The  boundary  conditions  were  taken  to  be  that  at  time 
zero,  the  gas  occupied  a  volume  filling  the  space  from  the  surface  out  to  a 


distance  equal  to  ^Tq/M) 


1/2 


times  the  pulse  length.  Here  M  is  the  mass  of 


the  particles  in  the  gas  cloud.  The  velocity  of  the  expanding  edge  at  time 
zero  was  taken  as  (kT^/M)^^.  and  at  T  =  O.  T^  -  MX  d^X 

(J  O - T 

3  k  dr 
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In  Figure  IE- 8  we  show  the  results  for  pulse  shapes  at  a  position  1  cm  from 
the  target  using  the  uniform  density  model.  These  results  were  obtained 
by  solving  equation  (III- 6)  numerically  with  the  given  boundary  conditions, 
and  give  the  pulse  shapes  due  to  the  arrival  of  mass  28  ions  at  a  position  1 
cm  from  the  target.  Qualitatively,  these  results  are  very  similar  to  those 
shown  in  Figure  V-  1  of  the  last  semi-  annual  report  where  we  presented 
analogous  results  obtained  cn  the  density  gradient  model.  The  pulse  shapes 
shown  in  Figure  ni-8  indicate  a  vertical  rise  and  very  sharp  peaks  for  the 
pulses  in  the  uniform  density  model.  Quantitatively,  both  the  uniform 
density  model  and  the  density  gradient  model  gives  pulses  that  rise  steeply 
and  decay  relatively  slowly.  In  addition,  both  models  give  approximately  the 
same  relative  change  as  the  assumed  value  of  T^changes.  As  Tq  increases, 
the  pulses  rise  and  decay  faster  and  show  decreased  widths.  The  results  for 
the  linear  density  gradient  model,  as  would  be  intuitively  expected,  give 
somewhat  more  rounded  peaks  which  seem  more  physically  reasonable, 
compared  to  the  very  sharp  peaks  for  the  pulses  obtained  on  the  uniform  density 
model.  Accordingly,  our  later  work  will  consider  only  the  density  gradient 
model  and  we  will  not  carry  through  the  analogous  derivation  for  the  uniform 
density  model. 

The  treatment  described  above  and  in  the  last  semi-annual  report  for  the 
free  expansion  model  gives  some  interesting  and  suggestive  results.  How¬ 
ever,  the  boundary  conditions  assumed  for  this  model  are  subject  to  some 
question.  For  example,  we  assume  that  the  thickness  of  the  sheet  of 
the  blowoff  material  at  time  zero  is  equal  o  (kT^  / M)  '  times  the  pulse 
length.  This  value  may  be  inexact  because  of  space  charge  effects  as  we  have 
seen  above  in  Section  III- A.  Also,  the  solution  derived  on  the  basis 
of  these  boundary  conditions  gives  an  asymptotic  value  for  the  velocity 
of  expansion  which  is  not  the  same  as  the  value  predicted  by  a  more  rigorous 
treatment  of  the  f  r  e  e  expansion.  The  details  of  the  asymptotic  solu¬ 
tions  for  the  .expansion  of  gas  into  a  vacuum  have  been  derived  by 

(17) 

other  workers,  for  example  E.  H.  Wedemeyer.  The  results  of  Wedemeyer 
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TIME  AFTER  LASER  PULSE  (MICROSECONDS  ) 

Figure  III- 8  Calculated  mass  28  particle  flux  at  r  position  1  cm  from  the  target 
based  on  free  expansion  model  with  uniform  density  assumption. 


can  be  applied  to  the  case  of  our  one -dimensional  geometry.  Tt  should  be 

noted  that  this  problem  is  well  known.  Wedemeyer's  results  agree  with  the 

results  obtained  by  others.  The  particular  thing  we  are  most  interested  in 

is  the  asymptotic  velocity  for  the  adiabatic  free  expansion  of  a  cloud  of  gas 

into  a  vacuum.  If  we  use  Wedemeyer's  result  we  obtain  an  asymptotic  value 

l/2 

of  the  expansion  velocity  equal  to  2  (  y  kT~/M)  /  (y  -  1).  For  y  =  1.  4, 
a  value  representative  oi  diatomic  gases  such  as  CO,  one  obtains  5.  91(kT^/M)  '  , 

and  for  y  =  1.67,  a  value  representative  of  monatomic  gases 

1/2 

such  as  alkali  metals,  one  obtains  3.  87  (kT~/M)  '  .  In  contrast,  on  the 

density  gradient  model  we  obtain  (7kT~/M)  ,  and  on  the  uniform  density 

l/2  ^ 

model  we  •  btain  ^kT^/M)  '  for  the  asymptotic  velocity.  Since  the  density 
gradient  model  gives  more  physically  realistic  pulse  shapes,  we  see  that  our 
previous  treatment  seems  to  indicate  asymptotic  expansion  velocities  that  are 
low  by  a  factor  of  2,  and  energies  for  the  directed  expansion  that  are  low  by 
a  factor  of  4. 

The  results  of  our  density  gradient  model  for  CO  and  the  more  correct  result 
given  by  Wedemeyer  agree  only  if  y  =  2.  09.  This  is  an  unreasonably  large 
value  of  y  for  CO.  Accordingly,  we  have  modified  our  treatment  to  preserve 
its  essential  features  but  to  give  an  accurate  value  for  the  expansion  velocity. 

These  two  difficulties  mentioned  above  can  both  be  dealt  with  by  using  a 
different  set  of  boundary  conditions.  Since  the  discussion  described  above  in 
Section  III-A  on  the  thickness  of  the  sheet  of  electrons  in  the  blowoff  material 
means  that  the  thickness  is  in  doubt,  we  discard  this  as  a  boundary  condition 
and  instead  use  the  asymptotic  value  of  the  velocity  of  the  expanding  edge. 
Some  results  obtained  by  numerical  integration  of  equation  (III- 6)  with  this 
boundary  condition  are  shown  in  Figures  III-9  and  III-10,  Figure  III- 9  shows 
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TIME  AFTER  LASER  PULSE  (MICROSECONDS) 

Figure  III-  9  Calculated  mass  28  particle  flux  1  cm  from  target  based  on 

asymptotic  value  for  free  expansion  with  density  gradient  assumption. 


300,00 


TIME  AFTER  LASER  PULSE  (MICROSECONDS) 

Calculated  mass  23  particle  flux  1  cm  from  target.  Based  on  asymptotic 
value  for  free  expansion  with  density  gradient  assumption. 


the  pulse  shape  obtained  at  a  collector  1  cm  distant  from  the  target  as  a  function 
of  time  for  the  free  expansion  of  a  gas  of  mass  28  molecules.  The  original 
assumed  temperature  Tq  is  a  parameter  in  these  curves.  These  results 
implicitly  assume  the  density  gradient  model.  This  figure  may  be  compared 

(4\ 

with  Figure  V-l  of  the  last  semi-annual  report.  *'  We  see  that  using  the 
asymptotic  velocity  of  expansion  as  a  boundary  condition  gives  curves  that 
generally  have  the  same  shape  as  those  obtained  with  the  earlier  formulation. 
However,  the  asymptotic  velocities  are  higher  and  the  corresponding  energies 
are  also  higher  for  a  fixed  initial  temperature.  The  material  begins  arriving 
at  a  collector  1  cm  from  the  target  at  earlier  time  and  also  the  pulses  are 
somewhat  sharper  with  shorter  width  at  half  height.  Figure  III -10  shows  similar 
results  for  pulse  shapes  for  mass  23  particles.  Here  because  of  a  larger 
value  of  y  »  the  corresponding  energies  are  lower  than  for  the  case  of  mass  28 
particles.  The  dominating  influence  here  is  the  (  y  ~  1 )  term  in  the  denominator 
of  the  asymptotic  velocity. 

For  mass  23  particles  the  energy  corresponding  to  the  velocity  of  the  expanding 
edge  can  correspond  to  a  particle  energy  as  high  as  190  eV  for  an  assumed  Tq 
value  of  300,  000°K.  If  we  had  to  explain  the  entire  energy  of  190  eV  on  a  purely 
thermal  basis,  it  would  require  temperatures  of  the  order  of  2  million  degrees  K. 
Thus  we  see  that  invoking  a  free  expansion  after  an  original  heating  does  reduce 
the  requirements  on  the  heating.  However,  even  heating  to  300,  000°K  requires 
additional  mechanisms  which  are  as  yet  unexplained.  We  saw  in  Section  III-A 
some  considerations  relating  to  the  use  of  the  Saha  equation  and  conditions  under 
which  we  may  obtain  neating  by  absorption  of  the  laser  light  in  the  inverse 
Bremsstrahlung  process,  These  results  indicate  possibly  a  higher  heating 

(2) 

rate  than  that  which  we  calculated  in  earlier  considerations  of  the  same  subject. 

The  results  of  the  present  calculation  on  the  free  expansion  model  using  the 
asymptotic  value  of  the  expansion  velocity  as  a  boundary  condition  are  summed 
up  in  Table  III -2  where  we  present  the  velocities  and  corresponding  particle 
energies  for  various  masses  and  original  temperatures  Tq.  We  see  from  this 
table  that  the  particle  energies  predicted  can  be  of  the  same  order  of  magnitude 
as  the  ion  energies  observed  in  the  time-of-flight  spectrometer. 
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TABLE  III -2 

ASYMPTOTIC  SOLUTION  FOR  DENSITY  GRADIENT  MODEL 
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In  summary,  the  picture  that  can  be  deduced  from  this  analysis  is  consistent 
with  the  high  energy  ion  pulses  being  produced  by  expansion  of  a  heated  gas 
cloud  provided  that  the  desorbed  gas  can  be  heated  rapidly  during  the  laser 
pulse.  This  leads  to  energies  of  expansion  with  superimposed  random 
components  of  the  same  order  of  magnitude  as  we  actually  observe  in  the 
time -of -flight  spectrometer,  and  also  to  pulse  shapes  similar  to  those  obtained 
experimentally.  We  shall  consider  some  of  the  speculations  concerning  a 
synthesis  of  the  processes  occurring  in  the  blowoff  material  as  derived  from 
both  the  time -of -flight  spectrometer  and  the  quadrupole  spectrometer  measure¬ 
ments  in  Section  III-E,  and  will  consider  the  mechanism  of  heating  that  could 
occur. 

Further  work  on  the  expansion  of  the  laser  produced  blowoff  material  will 
include  the  following  steps: 

1.  Inclusion  of  the  effect  of  an  accelerating  potential.  In  our  actual 
measurements  on  the  time -of -flight  spectrometer,  ions  are  accelerated 
across  the  gap  by  an  applied  field  of  the  order  of  kilovolts.  This 
applied  field  also  repels  the  electrons  into  the  target  and  affects  the 
space  charge  sheath  near  the  target.  It  will  be  necessary  to  include 
the  effect  of  the  space  charge  and  also  the  motion  of  the  entire  cloud 

of  ions  away  from  the  target  surface.  In  this  formulation  the  density 
gradient  model  as  we  have  set  it  up  previously  will  no  longer  be 
applicable. 

2.  Inclusion  of  drift  time  down  the  drift  tube  in  the  time -of -flight  spectro¬ 
meter.  When  this  drift  time  is  included,  we  will  obtain  pulse  shapes 
that  can  be  compared  directly  to  experimental  pulse  shapes. 
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3.  Inclusion  of  possible  absorption  in  the  plasma.  Equation  (IiI-8)  can 
be  modified  to  include  a  term  for  possible  absorption  in  the  plasma. 

The  right  hand  side  will  be  non-zero.  The  conjectures  described 
previously  about  the  heating  rate  may  be  higher  than  we  considei  jd 
it  to  be  in  earlier  calculations  and  will  give  us  insight  into  how  this 
whole  routine  can  be  modified  to  include  the  effect  of  heating  during 
the  early  stages  of  the  expansion.  This  whole  formulation  would  then 
have  to  be  solved  numerically  in  order  to  give  a  complete  picture  of 
how  the  blowoff  material  heated  and  expanded. 

4.  In  the  computer  runs  in  which  the  results,  such  as  are  shown  in 
Figures  III-9  and  III-10,  were  calculated  we  found  that  the  velocity  of 
the  expanding  edge  very  rapidly  approached  the  asymptotic  velocity. 

In  fact,  within  only  a  few  nanoseconds  after  the  material  began 
expanding  it  had  reached  the  asymptotic  velocity  and  most  of  the 
distance  from  the  target  to  the  distance  1  cm  away  was  covered  at 

this  asymptotic  velocity.  This  means  that  the  picture  of  a  30  nanosecond 
long  pulse  in  which  the  blowoff  material  is  produced,  and  raised  to  a 
temperature  Tq,  and  then  expands,  cannot  be  realistic  since  the  time  of 
transition  from  the  original  state  to  the  asymptotic  state  is  shorter  than 
the  time  we  consider  the  material  as  remaining  in  its  initial  state.  Our 
future  work  will  take  this  fact  into  account  and  consider  simultaneous  heat¬ 
ing  and  expansion. 

D.  ELECTRON  EMISSION 

Calculations  have  been  carried  out  on  the  electron  emission  from  a  laser 

heated  surface.  These  calculations  are  essentially  a  return  to  work  done 
(18) 

some  time  ago.  v  We  had  found  a  discrepancy  between  calculated 

temperatures  based  on  thermodynamic  considerations  and  observed 

temperatures  as  deduced  from  the  electron  emission  in  the  laser-heated 
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spot  using  Richardson's  equation.  In  particular,  as  seen  from  Figures  3 

(18) 

and  4  of  reference  '  ,  the  experimentally  deduced  temperatures 

consistently  ran  several  hundred  degrees  higher  than  the  theoretically 
calculated  temperatures.  At  the  time  we  hypothesized  that  the  reason  for 
this  discrepancy  was  that  the  power  density  in  a  laser  beam  was  not  uniform. 
This  non -uniformity  gave  local  hot  spots  in  the  laser  beam  and  contributed 
largely  to  the  emission  because  of  the  non-linear  dependence  of  the  electron 
emission  on  temperature.  These  local  hot  spots  dominate  the  emission  to 
the  extent  that  it  makes  the  average  temperature  look  higher  than  it  actually 
is.  Most  of  the  electron  emission  arises  from  the  area  near  the  center  of 
the  laser  beam  where  the  temperature  is  highest.  Even  though  this  is  a 
relatively  small  fraction  of  the  area,  the  dependence  of  the  electron  emission 
is  so  strongly  a  function  of  temperature  that  when  the  total  electron  emission 
is  divided  by  the  total  area,  one  obtains  a  value  for  the  average  emission 
that  corresponds  to  a  temperature  higher  than  the  average  temperature  of 
the  surface. 

In  the  present  work  we  used  the  spatial  distribution  of  power  in  the  laser  pulse 
as  described  in  a  previous  report,  ^  and  set  up  a  calculation  for  a  laser  power 
distribution  with  this  shape.  The  treatment  here  is  similar  to  the  calculation 
of  the  desorption  of  carbon  monoxide  from  a  tungsten  surface  as  described  in 

/  n  \ 

a  previous  report.  '  In  that  work  we  started  with  a  known  rate  of  desorption 
of  carbon  monoxide  as  a  function  of  temperature.  The  temperature  was  then 
calculated  as  a  function  of  position  and  time  for  a  laser  pulse  with  specified 
spatial  and  temporal  shapes,  and  numerically  integrated  to  obtain  the  total 
desorption  of  carbon  monoxide.  In  the  calculation  it  was  assumed  that  there 
was  essentially  no  transverse  flow  of  heat  in  the  laser-heated  spot.  In  other 
words,  the  transverse  dimensions  of  the  spot  were  much  greater  than  the  depth 
of  penetration  of  the  heat  during  the  time  of  the  laser  pulse.  This  condition 
is  reasonably  well  approximated  for  our  work.  It  means  that  the  rise  in 
temperature  at  a  position  in  the  laser  heated  spot  is  directly  proportional  to 
the  input  power  density  at  that  position. 
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For  the  electron  emission  as  a  function  of  temperature  at  a  particular 
position  in  the  laser  spot  we  used  Richardson's  equation.  For  the  laser 
power  density  as  a  function  of  position  we  used  the  distribution  described 
previously.  ^  As  a  function  of  the  radius  r  from  the  center  of  the  spot 
and  the  polar  angle  9, 

P  (r,  0)  =  P  [  1  -  - - -  ] 

u  a  (n-sin  26) 


Pq  is  the  power  per  unit  area  at  the  center  of  the  illuminated  area. 


Here  the  approximate  values  of  n  and  ot  are  n  =  4,  and  =  0.  075  cm.  The 

variables  in  the  temperature  T  (r,  9,  t)  may  be  separated  according  to  the 
formulation  described  above  in  which  the  temperature  rise  at  a  particular 
position  is  linearly  proportional  to  the  power  density  at  that  particular 
position  under  the  assumption  that  there  is  no  transverse  flow  of  heat.  In 
other  words,  T(r,  9,  t)  equals  P{r,  9)  times  the  temperature  as  a  function  of 
time  for  a  unit  power  uniform  spatial  distribution  laser  pulse  with  the 
desired  temporal  shape.  With  this  formulation  the  electron  current  as  a 
function  of  time  is: 
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In  this  equation  A  is  the  constant  in  Richardson's  equation,  T^(t)  is  the 
temperature  for  a  spatially  uniform  laser  pulse  of  unit  power  with  the 
desired  temporal  shape,  Pq  is  the  peak  laser  power  per  unit  area  at  the 
center  of  the  laser  spot  at  the  time  at  which  the  laser  power  is  highest, 

0  is  the  work  function  of  the  surface,  k  is  Boltzman's  constant,  Tq  is  the 
initial  temperature,  and  the  other  variables  are  defined  above.  This 
equation,  by  suitable  manipulation,  can  be  reduced  to  the  following  form: 

J(t)  =  APq  2Tl  2  (t)  a  *  ir  (2n2  +  1) 


+  (  l+^> 


+  ?2(t?e  +^}  E  (a(t))  ] 

6  6  6  1 


Here  a(t)  =  0/k  (Tq  +  PqT^O),  and  E^aft))  is  the  exponential  integral. 

(15) 

We  have  used  values  for  T^(t)  taken  from  earlier  calculations,  for  a  laser 

pulse  shape  actually  representative  of  that  which  was  used  in  the  electron  emission 

(18) 

experimental  work.  v  The  above  equation  was  then  evaluated  on  a  computer 
using  a  numerical  routine  for  evaluation  of  (a(t))  with  the  values  of  T^(t) 
stored  in  memory.  Some  results  are  shown  in  Figure  III -11,  where  we  see  the 
electron  emission  from  tungsten  as  calculated  on  this  basis  of  nonuniform 
spatial  power  distribution.  The  parameter  is  the  laser  power  at  the  peak  value 
in  time  divided  by  the  area  of  the  illuminated  spot. 


One  may  take  results  such  as  shown  in  Figure  III- 11  and  use  the  known  values 
for  thermionic  emission  as  a  function  of  temperature  and  deduce  an  apparent 
temperature  corresponding  to  the  average  value  of  electron  emission  over  the 
whole  spot.  This  is  the  quantity  that  we  compare  to  an  experimentally  deduced 
temperature  from  the  laser  emission.  The  results  are  shown  in  Figure  III  - 1 2 
for  t’noriated  tungsten.  This  is  essentially  a  modification  of  Figure  3  of  the 
reference  (18).  The  present  version  of  this  figure  gives  a  postulated  laser 
power  as  a  function  of  time.  This  is  the  total  laser  power  divided  by  the  total 
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TOTAL  ELECTRON  EMISSION  DIVIDED  BY  ILLUMINATED  AREA  (AMP/CM*) 
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TIME  AFTER  START  OF  LASER  PULSE  ( NANOSECONDS  l 

Figure  III- 1 2  Apparent  temperature  of  thoriated  tungsten  surface  as  derived  from  calculation 
and  from  experiments  on  electron  emission  for  a  laser  pulse  of  the  given  shape. 
Right  scale  is  total  laser  power  divided  by  illuminated  area. 


area  of  the  spot.  Figure  III-12  shows  the  experimental  temperature  as 

derived  from  the  measured  values  of  the  electron  emission,  and  also  two 

theoretical  values,  first  the  theoretical  value  obtained  from  assuming  that 

the  laser  power  is  uniform  over  the  spot,  and  second  the  theoretical  value 

as  modified  for  the  nonuniform  temperature  distribution  as  described  above. 

We  see  that  in  this  particular  power  range  (peak  power  around  12  megawatts/ 

2 

cm  average  value)  using  the  spatial  distribution  actually  representative  of 
our  laser  pulse  gives  a  temperature  value  several  hundred  degrees  higher 
than  the  temperature  calculated  assuming  a  uniform  spatial  distribution. 

This  calculated  temperature  is  in  quite  good  agreement  with  the  experi¬ 
mentally  determined  temperature  as  a  function  of  time. 

For  laser  powers  of  approximately  10  to  15  MW/cm  good  agreement  is 
obtained  between  the  presently  modified  theory  and  the  experiment.  However, 
at  higher  powers  the  theoretical  values  are  rising  faster  than  the  experimental 
values.  This  occurs  because  the  center  of  the  spot  becomes  very  hot,  and 
because  of  the  non-linear  dependence  of  electron  emission  cn  temperature. 
Thus  the  values  calculated  on  the  present  theoretical  formulation  become 
much  higher  than  the  experimentally  determined  values  at  laser  powers  of  the 
order  of  20  to  25  MW/cm  .  This  can  be  explained  by  saturation  of  the  experi¬ 
mentally  obtained  current  as  the  current  becomes  very  high.  Space  charge 

effects  limit  the  experimental  current  as  the  temperature  at  the  center  of  the 

(18) 

spot  rises.  In  the  earlier  work'  we  had  a  collecting  voltage  high 

°  2 

enough  to  overcome  space  charge  effects  at  laser  powers  around  15  MW/cm  . 
However,  at  the  higher  laser  powers  at  the  center  of  the  spot  we  can  have 
enough  space  charge  effect  so  that  we  no  longer  can  pull  all  the  electrons 
over  to  the  collector. 
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In  summary,  using  experimentally  determined  spatial  distributions  for 

the  laser  power  density  gives  good  agreement  between  theoretical  and 

experimentally  determined  values  of  the  temperature  at  the  low  end  of  the 

range  of  laser  powers  used  in  the  electron  emission  experiments.  However, 

2 

at  higher  powers  of  the  order  of  20  to  25  MW  /  cm  ,  at  the  upper  end  of 

the  range  of  laser  power  densities  used  in  the  experimental  electron  emission, 

there  still  remain  discrepancies  between  theory  and  experiment. 

E.  PARTICLE  EMISSION  IN  LASER- SURFACE  INTERACTIONS 

A  coherent  picture  of  the  entire  system  of  particles  produced  in  a  laser- 
surface  interaction  as  obtained  from  the  combined  results  of  measurements 
on  the  time-of-flight  spectrometer  and  the  quadrupole  spectrometer  may 
be  synthesized  to  explain  the  range  of  phenomena  observed. 

The  ideas  presented  in  this  subsection  are  preliminary  and  speculative. 

We  do  not  have  proof  of  all  the  hypotheses  given  here.  They  do  appear 
to  be  physically  appealing  and  are  consistent  with  our  experimental  results, 
but  further  work  will  be  required  to  provide  more  definite  proof. 

As  the  surface  is  heated,  we  get  emission  of  neutral  gases  such  as  CO,  of 
alkali  metals  emitted  as  ions,  for  example,  sodium  and  potassium  ions,  and 
also,  at  our  power  densities,  a  relatively  small  amount  of  the  base  material 
such  as  tungsten.  This  is  accompanied  by  a  fairly  large  amount  of  thermionic 
electron  emission.  The  electrons  set  up  a  space  charge  sheath  near  the 
surface  in  a  fairly  short  time.  The  equilibrium  time  during  which  this 
sheath  is  set  up  can  be  very  short  because  of  the  high  mobility  of  the  electrons, 
and  we  can  reach  a  steady  state  even  during  the  time  of  the  laser  pulse.  This 
local  region  of  high  charge  density  can  lead  to  high  rates  of  heating. 
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In  particular,  we  look  on  the  sodium  ions  as  being  rapidly  heated  by 
absorption  of  the  laser  light  in  the  inverse  Bremstrahling  process  in  this 
localized  region  of  high  charge  density.  Most  of  the  material  does  not 
heat  to  high  temperatures.  This  is  evidenced  by  the  fact  that  the  CO  is 
emitted  mainly  as  undissociated  molecules. 

We  consider  the  main  mechanism  of  energy  transfer  to  be  through  electron 
collisions.  We  consider  the  electron  collisions  to  be  rapid  enough  that 
thermal  equilibrium  can  be  established  on  a  local  basis  between  the  electrons 
and  the  ions,  but  the  interaction  of  electrons  and  neutral  particles  to  be  slow 
enough  that  essentially  little  heating  of  the  neutral  molecules  can  occur. 
Thereafter  we  can  have  elastic  collisions  between  energetic  sodium  ions 
and  the  thermal  CO  molecules.  This  will  lead  to  some  dissociation  of  the 
CO  but  there  will  also  be  elastic  collisions  in  which  high  energy  CO  molecules 
will  be  produced. 

Relative  numbers  of  the  various  components  will  have  to  agree  with  this 
picture.  In  particular,  we  should  see  much  more  thermal  CO  than  any  other 
component.  Also,  collisions  with  the  hot  sodium  ions  should  produce 
dissociation  of  some  of  the  CO.  Thus,  we  should  see  some  low  energy 
carbon  and  oxygen  produced  in  the  dissociation  interactions.  We  should  also 
see  a  smaller  amount  of  high  energy  CO  molecules.  Experimentally  the 
numbers  of  high  energy  CO  molecules  that  we  observe  is  somewhat  in  doubt 
because  of  uncertainty  about  the  secondary  emission  coefficient  of  the  photo¬ 
cathode  for  energetic  molecules.  However,  the  results  are  not  inconsistent 
with  their  being  many  fewer  high  energy  CO  molecules  than  thermal  CO 
molecules.  This  is  the  case  if  the  secondary  emission  coefficient  approaches 
unity.  This  picture  is  also  supported  by  the  observation  that  the  high  energy 
CO  typically  has  about  half  as  much  energy  as  the  high  energy  ions,  as  might 
be  expected  if  they  were  produced  in  elastic  collisions  with  energetic  alkali  ions. 
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This  picture  synthesizes  the  main  features  of  the  observations  made  on  two 
different  types  of  spectrometers.  It  provides  a  plausible  heating  mechanism 
for  the  energetic  ions  that  are  in  fact  observed,  but  allows  most  the  material 
to  remain  cool.  The  mechanism  for  the  high  energy  neutral  molecules 
observed  is  assumed  to  be  elastic  scattering  from  high  energy  alkali  ions. 

The  directed  energy  of  the  center  of  mass  of  the  entire  system  is  assumed  to 
come  from  an  adiabatic  free  expansion. 

This  picture  must  be  regarded  as  preliminary,  but  it  does  combine  enough 
features  to  explain  the  results  obtained  on  both  types  of  spectrometers,  and 
has  physically  appealing  reasons  to  back  it.  The  great  question  is  whether 
there  is  enough  heating  possible  to  heat  the  ions  to  the  required  high  tempera¬ 
tures.  If  so,  we  can,  as  we  saw  in  Section  III-C,  explain  the  pulse  shapes, 
the  directed  energies,  and  the  superimposed  random  components  of  energy 
for  the  alkali  ion  pulses  observed  in  the  time -of -flight  spectrometer.  The 
reasoning  given  in  Section  III -A  will  be  developed  further  to  obtain  quantitative 
estimates  of  the  heating  for  our  experimental  conditions  in  order  to  answer 
this  question. 


-  81  - 


SECTION  IV 


SUMMARY 

The  work  described  in  this  report  deals  with  a  number  of  related  topics 
relevant  to  the  laser-surface  interaction.  The  experimental  and  theore¬ 
tical  work  complement  each  other  and  will  lead  to  an  understanding  of 
the  nature  of  the  processes  occurring  when  a  high  power  laser  begun 
strikes  a  metgillic  surface.  Some  of  the  topics  described  here  are  ex¬ 
tensions  of  measurements  or  calculations  which  have  been  carried  out 
previously  as  part  of  this  continuing  effort,  for  exaimple,  measurements 
on  the  ion  and  neutrgil  molecule  emission  from  sodium  targets,  and  cal¬ 
culations  on  the  particle  pulse  shapes  resulting  from  a  free  expansion  of 
the  heated  blowoff  material.  Other  topics  are  introduced  in  this  report 
for  the  first  time;  for  example,  measurements  of  the  angular  distribu¬ 
tion  of  the  ionic  emission,  and  calculation  of  thermomechanical  stress 
profiles  in  the  target. 

At  this  point  we  stand  at  a  trguisition  phase  in  this  work,  with  earlier  types 
of  measurements  essentially  being  completed,  and  new  types  of  measure¬ 
ments  yielding  preliminary  results.  The  earlier  work,  on  ion  and  neutral 
gas  emission  from  various  surfaces  at  relatively  low  laser  powers  (of 
the  order  of  50  MW /cm  )  ,  which  relied  mainly  on  the  quadrupole  and 
time -of- flight  mass  spectrometers,  is  being  replaced  by  other  types  of 
measurements  on  new  equipment  designed  to  elucidate  specific  features 
of  the  laser- surface  interaction,  for  example,  the  angular  distribution 
of  ion  emission.  The  capability  of  extending  the  megLSurements  tc.  a  high¬ 
er  range  of  laser  powers  has  also  been  developed.  We  have  also  reached 
the  stage  where  the  various  results,  both  theoretical  and  experimental, 
can  be  integrated  to  provide  a  complete  picture  of  the  phenomena  occur- 
ing  in  the  laiser- produced  blowoff  material,  whereas  the  earlier  work 
tended  to  consider  one  specific  phenomenon  at  a  time. 
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The  specific  results  described  in  this  report  include  the  following  as 
the  most  significant: 

1.  Preliminary  indications  of  an  anisotropy  in  the  ion  emission, 
with  most  ions  emitted  in  the  direction  normal  to  the  surface. 

2.  Measurements  on  sodium  targets,  indicating  large  sodium 
ion  and  hydrocarbon  fragment  emission. 

3.  Refinement  of  calculations  on  the  free  expansion  of  the  blow- 
off  material,  indicating  that  the  results  for  an  initial  tempera- 

5o 

ture  around  10  K  are  similar  to  the  pulse  shapes,  including 
the  directed  velocity  component,  of  the  observed  ion  emission. 

4.  Investigations  of  the  space  charge  effects  and  ionization  levels 
in  the  blowoff  material,  with  implications  relative  to  the 
possible  heating  of  this  material  by  the  laser  light. 

5.  Preliminary  formulation  of  an  integrated  picture  to  describe 
experimental  results  from  both  the  time-of- flight  and  quad¬ 
ruple  spectrometers,  in  which  heating  occurs  in  a  plasma 
of  alkali  ions,  while  the  neutral  gas  remains  cool. 

6.  Calculation  of  thermomechanical  stress  profiles  in  the  target 

material  due  to  absorption  of  the  laser  light,  in  which,  for  our 

experimental  conditions,  the  main  feature  is  a  sharp  pressure 

8  2 

spike  of  the  order  of  6  x  10  dyne/cm  propagating  into  the 
target. 

For  the  future,  plans  call  for  completing  and  extending  the  work  describ¬ 
ed  above.  Specifically  we  shall  complete  the  investigation  of  the  angular 
distribution  of  the  ion  emission  and  carry  out  a  critical  experiment  to 
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determine  the  rolo  of  the  absorption  of  the  laser  light  in  heating  the 
blowoff  material.  We  shall  extend  the  measurements  into  a  higher  range 
of  laser  powers  than  has  been  used  previously.  The  plasma  heating 
mechanism  will  be  investigated  in  the  light  of  the  considerations  of  ioni¬ 
zation  levels  in  the  blowoff  material  as  described  above.  These  results 
will  be  integrated  with  the  work  on  expansion  of  the  blowoff  material 
in  order  to  analyze  the  formulation  for  the  picture  of  the  particle  emission 
as  described  in  Section  III-E.  The  calculations  on  the  shock  waves  gen¬ 
erated  in  the  targets  will  be  extended  to  include  other  target  materials. 
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